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'Satellite  observations  have  the  ability  to  determine  meteorological 
conditions  over  inaccessible  or  unfriendly  territory.  The  Defense 
Meteorological  Satellite  Program  (DMSP)  has  two  satellites  (Bird 
8531  and  9532)  that  provide  infrared  radiance  data  on  a routine 
basis  for  operational  applications  at  the  Air  Force  Global  Weather 
Center  (AFGWC).  Multi-spectral  infrared  radiometers  operating  at 
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identical  wavelengths  on  each  of  the  satellites  are  used  to  interro- 
gate different  depths  in  the  atmosphere  to  derive  air  temperature 
profile  information.  AFGWC  has  a sophisticated  computer  program 
that  operates  upon  satellite  radiance  observations  using  their  12- 
hour  forecast  temperature  fields  as  a first  guess  to  derive  a 
measure  of  the  temperature  profile.  Errors  are  especially  noticeable 
in  the  lower  troposphere  near  the  earth's  surface  and  in  the  tropo- 
pause  region.  This  report  describes  an  analytical  and  a statistical 
approach  to  isolate  the  largest  error  sources  and  to  suggest  promising 


avenues  for  operation  improvement.  A research  version  of  SOUNDER, 
the  AFGWC  operational  computer  prjSgram,  was  formulated  and  used 
to  systematically  analyze  physical,  cbijiputational , and  hardware 
contributions  to  errors  in  the  retrieval  oXtemperature  profiles.  It 
was  found  that  radiance  observations  from  Bird  8531  were  higher 
(warmer)  than  for  Bird  9532  in  seven  of  the  eight  sensor  channels. 
The  present  operational  criteria  used  to  eliminate  low  level  clouds 
or  hot  surface  are  so  restrictive  that  many  of  the  most  significant 
types  of  atmospheric  temperature  profiles  are  also  eliminated. 

Water  vapor  contributions  to  errors  can  be  large  in  the  tropics, 
worthy  of  consideration  in  mid-latitudes,  and  insignificant  in  the 
polar  regions.  Temperature  retrieval  errors  due  to  water  vapor  are 
most  pronounced  at  the  7 00  mb  level.  Underlying  surface  effects 
are  not  properly  taken  into  account.  These  surface  effects  provide 
the  greatest  single  error  source  for  AFGWC 's  operational  routines 
for  satellite  temperature  retrievals  in  the  lower  troposphere.  It  is 
shown  that  these  surface  effects  have  an  influence  on  all  sensor 
channels  and  are  related  to  the  low  level  atmospheric  stability 
conditions.  Statistical  models  were  developed  and  compared  with 
the  minimum  information  technique  in  current  use  at  AFGWC  for 
temperature  retrievals.  Significant  error  reductions  are  possible 
using  the  statistical  methods.  It  appears  that  statistical  coefficients 
derived  from  one  week  of  data  can  be  used  with  confidence  during 
the  subsequent  week  to  derive  temperature  profiles  from  multi- 
channel satellite  radiance  data. 


Qualified  requestors  may  obtain  additional  copies  from  the  Defense 
Documentation  Center.  All  others  should  apply  to  the  National 
Technical  Information  Service. 
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1 INTRODUCTION 


1.1  GENERAL 

1.1.1  Overview 

Meteorological  observing  station  density  increases  have  improved 
weather  analysis  and  forecast  accuracy.  Many  uninhabited  regions  of  the 
world  exist,  especially  over  ocean  and  polar  regions,  that  can  be  in- 
terrogated remotely  from  satellites  to  obtain  weather  information.  The 
Defense  Meteorological  Satellite  Program  (DMSP)  has  provided  routine 
observations  of  weather  phenomena  on  a scale  never  before  attainable. 
Improvements  in  the  aviation  weather  variables  of  cloud,  fog,  and 
significant  weather  depiction  have  been  obtained.  The  advent  of  multi- 
frequency infrared  radiometric  sensors,  operating  on  absorption/emission 
bands  of  atmospheric  constituents,  has  opened  new  horizons  for  remotely 
probing  the  vertical  structure  of  the  atmosphere.  This  has  been  especially 
true  for  remotely  probing  the  vertical  temperature  profile  which  plays  an 
important  role  in  specifying  atmospheric  stability  that  is  so  important  to 
numerical  weather  prediction.  Performance  of  present  methods  for  remotely 
deriving  temperature  profiles  from  satellites  has  been  questioned.  Even 
for  cloud  free  lines  of  sight,  large  errors  exist  in  attempting  to  retrieve 
air  temperatures  from  satellite  radiance  data.  These  large  errors  are 
especially  noticeable  in  the  lower  troposphere  near  the  surface  and  in 
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the  tropopausc  region  separating  the  troposphere  from  the  stratosphere. 

It  is  the  source  of  these  large  errors  and  the  identification  of  promising 
corrective  methods  that  received  emphasis  in  this  study. 

Two  methods  of  approach  taken  in  this  study  included  analytical 
and  statistical  techniques.  An  analysis  was  made  of  Air  Force  Global 
Weather  Central's  (AFGWC)  "SOUNDER"  computer  program  package  for 
operationally  retrieving  temperature  profiles  from  12-hour  forecast  and 
satellite  data  inputs.  A Geo-Atmospherics  Corporation  (GAG)  research 
version  of  SOUNDER  was  developed  and  used  to  perform  a sensitivity 
analysis  of  the  response  of  SOUNDER  to  accurate  and  contaminated 
data,  error  source  discrimination  criteria,  solution  convergence  pro- 
cedures, sensor  noise  and  errors,  water  vapor,  and  underlying  surface 
effects.  A statistical  analysis  was  made  on  a large  two  week  February 
197  5 northern  hemisphere  data  base  at  Air  Force  Geophysics  Laboratory 
(AFGL).  Data  files  of  cloud  free  satellite  soundings  and  "nearby  ground 
truth"  grid  point  temperature  analyses  were  generated  and  used  to 
derive  and  verify  stepwise  regression  equations  for  retrieving  air  tempera- 
tures from  satellite  data. 

Serious  sources  of  large  errors  were  identified,  especially  those 
due  to  underlying  surface  effects,  and  promising  methods  for  error  re- 
duction were  demonstrated.  Significant  reductions  in  the  rms  error  at 
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all  levels  in  the  atmosphere  were  obtained  with  the  statistical  models. 
The  remainder  of  this  report  provides  background  information,  details 
I of  SOUNDER,  discussion  of  analytical  and  satistical  results,  and  a 

I 
! 

list  of  conclusions  and  recommendations. 

1.1.2  Background 

The  problem  of  obtaining  vertical  temperature  profiles  from 
passive  satellite  sensors  has  been  the  object  of  much  study  for  almost 
20  years.  Kaplan  ^ ^ first  suggested  that  spectral  measurements  in 
the  IS^i  m CO2  band  could  be  used  for  temperature  soundings.  The 
first  experimental  verification  came  10  years  later  with  the  launch  of 
the  Nimbus  3 satellite.  Since  that  time,  a number  of  other  satellites 
capable  of  such  measurements  have  been  launched,  including  the  SSE 
package  utilized  in  the  DMSP.  ‘■ 

The  ability  to  perform  accurate  soundings  on  a world-wide 
basis  is  an  imporcant  element  in  meeting  the  desire  of  the  Air  Force 
to  achieve  a continuous  world-wide  weather  analysis  and  prediction 
capability.  Temperature  profiles  are  also  an  element  in  predicting 
likelihoods  of  cloud  type  and  cover. 
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1.1.3  Approaches  to  Temperature  Inversion 

Two  basic  approaches  are  used  for  temperature  inversion. 

These  may  be  classified  as: 

(1)  Atomistic  Approach 

(2)  Holistic  Approach 

The  atomistic  approach  is  based  on  the  use  of  micromodels. 

All  known  physics  is  used,  including  equations  of  radiative  transfer, 
hydrostatic  equation  and  Planck's  equation.  Mixing  ratios  of  the  various 
atmospheric  constituents  are  computed  and  explicit  calculations  are  made 
to  include  the  effects  of  temperature  on  liquid  water  content.  Computations 
are  based  either  on  physical  laws  or  empirically  derived  relations.  Modi- 
fications are  made  on  the  microsystem  level  in  order  to  obtain  a better 

( 2 ) 

microsystem  (subsystem)  model.  The  minimum-information  method  is  an 

example  of  an  atomistic  method  of  temperature  retrieval.  Other  examples 
are  Chahine's  method^  ^ ^ and  Smith's  nonlinear  iterative  method.^  ^ ^ 

The  holistic  approach  utilizes  an  information  model.  Models 
are  built  up  which  best  fit  the  needs  of  the  application  at  hand.  For 
the  temperature  inversion  problem,  holistic  methods  are  statistical  in 
nature  and  are  based  on  achieving  the  best  inversions  as  determined  via 
testing  on  an  independent  data  base.  Examples  of  holistic  inversion  methods 
are  the  statistical  eigenvector  method  (Smith  and  Woolf) ^ ^ ^ and  regression 
analysis. 

It  should  be  mentioned  that  some  methods  are  a combination  of 
atomistic  and  holistic  methods.  For  example,  the  "full  statistical" 
method  (Fritz  et  al.)^^  ^ uses  botli  transmittance  functions  and  temperature 
and  radiance  covariance  matrices.  It  should  be  pointed  out,  however,  that 
the  minimum-information  method  is  not  a holistic  approach  even  though  it  is 
usually  thought  of  as  a special  case  of  the  "full  statistical"  approach. 

This  is  due  to  the  fact  that  the  scalar  "noise-to-signal  power"  ratio  is 
selected  on  the  basis  of  expected  calibration  errors  and  expected  errors  in 
the  guessed  radiances  (Smith,  Woolf  and  Fleming).^  ^ ^ In  the  holistic 
approach,  this  parameter  would  be  chosen,  for  example,  to  minimize  the  mean- 
squared  prediction  errors  In  the  vertical  temperature  profiles,  using  an 
independent  data  base  for  performance  testing. 
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Based  on  this  discussion,  it  is  clear  that  the  essential 
difference  between  the  two  methods  is  the  way  in  which  models  are  built 
up.  Every  temperature  inversion  technique  requires  a model  which  relates 
the  measured  radiances  and  the  vertical  temperature  profile.  The  model 
may  or  may  not  include  other  parameters,  such  as  season,  latitude,  scanning 
pattern,  initial  forecast  field,  and  so  forth.  The  essential  difference 
in  the  methods  is  that,  in  the  atomistic  approach  the  model  is  developed 
based  on  physical  models  using  little  or  no  actual  data.  On  the  other  hand, 
the  holistically-derived  model  is  based  on  analysis  of  a significant  amount 
of  actual  data  for  which  independent  verification  may  be  obtained. 

Neither  the  atomistic  nor  the  holistic  approach  offers  a 
panacea  for  all  problems  and  the  choice  of  method  depends  upon  both  the 
nature  of  the  problem  and  the  availability  of  a representative  data  base. 

For  initial  development,  an  atomistic  approach  is  almost  mandatory,  since 
little  or  no  actual  data  is  available.  However,  as  data  becomes  available, 
the  use  of  a holistic  method  should  be  considered. 

The  development  of  methods  for  temperature  retrieval  has,  in 
fact,  followed  this  path.  The  most  widely  used  technique  is  the  minimum- 

information  technique,  used  in  conjunction  with  the  empirical  correction 

( 7 ) 

technique  of  Weinreb  and  Fleming.  Empirical  and  physical  models  are 

used  for  computation  of  the  transmlttances . Several  hybrid  techniques  were 

( R ) 

then  tested.  Fleming  and  Smith  compared  the  iterated  "full  statistical" 

method  to  the  iterated  minimum-information  method,  the  iterative  nonlinear 
( 9 ) 

method  of  Chahine  and  a general  nonlinear  iterative  method.  Experimental 
evaluation  indicated  that  the  "full  statistical"  method  gave  the  best  overall 
performance.  This  is  an  interesting  result  since  this  is  the  only  method, 
of  the  four  analyzed,  which  included  statistical  information  based  on  actual 
data  in  the  model. 

In  another  related  study,  Crosby  and  Weinreb^  ^compared  the 
performance  of  the  "full  statistical"  method  and  the  minimum-information 
method  for  December-January  data.  Their  results  indicated  that  minimum 
information  retrievals  were  uniformly  worse,  with  rms  differences  of  up  to 
1°K  in  the  region  between  200  and  300  mb  and  about  1/2°K  above  300  mb. 
However,  there  was  also  a degradation  of  about  1°K  at  800  mb. 
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( 12  ) 

Malkevich,  et  al.  used  empirical  orthogonal  functions  in 

a linear  regression  model  to  estimate  temperature  profiles.  The  empirical 
functions  were  found  by  using  the  kernel  function  of  the  integrated  form  of 
the  equation  of  radiative  transfer. 

Recently  there  has  been  an  increasing  interest  in  the  use  of 

statistical  methods.  Smith  and  Woolf ^ ^ ^used  a statistical  eigenvector 

method  for  estimation  of  vertical  profiles  of  temperature,  moisture  and 

cloudiness.  Comparisons  with  the  minimum- information  method  indicated  that 

the  rms  error  in  predicting  cloud-free  brightness  was  significantly  reduced 

( 12  1 

with  the  statistical  eigenvector  method.  Klein,  Kyle  and  Smith 
compared  the  minimum- inf ormation  and  statistical  eigenvector  methods  on 
operational  SSE  data  and  concluded  that,  on  the  basis  of  the  data  studied, 
neither  could  be  strongly  recommended  over  the  other.  The  statistical 
eigenvector  solution  was  better  for  tropical  soundings  and  in  the  lower 
troposphere.  At  mid-latitudes,  the  minimum-information  solution  was  better 

between  250  and  500  mb  but  poorer  elsewhere. 

(11) 

In  yet  another  study.  Smith  investigated  the  effect  of  mea- 

surement noise  on  both  the  minimum-information  and  statistical  eigenvector 
retrievals.  Both  were  found  to  be  very  sensitive  to  noise,  if  the  model  did 
not  account  for  the  noise.  The  sensitivity  v/as  substantially  reduced  when 
the  noise  was  accounted  for,  however. 


1.2  METHODS  OF  ANALYSIS 

1.2.1  Theoretical  Model  for  Temperature  Retrieval 
through  a Cloudless  Atmosphere 
For  a non-scattering  cloud-free  atmosphere  in  local  thermodynamic 
equilibrium,  the  integral  form  of  the  radiative  transfer  equation  is 


N(V)  = B[v,T(p  )]  T (v,p  ) 
s s 

- ' B|v.t(p)1 


(1.1) 
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where 

N(v)  = outgoing  spectral  radiance  within  a spectral  window 
centered  at  frequency  v 

B[v,T]  = Planck  radiance  for  temperature  T at  frequency  v 
T(v,p)  = transmittance  of  the  atmosphere  above  pressure  p 
x(p)  = arbitrary  monotonic  function  of  p 

p^  = atmospheric  pressure  at  the  earth's  surface 
The  objective  is  to  solve  (1.1)  for  T(p)  given  a finite  set 

of  radiance  measurements  N(v^),  NCv^), significant 

problem  in  (1.1)  is  the  nonlinear  dependence  of  B[v,T(p)]  on  T(p): 


B[v,T(p); 


exp 


s ] 

T(p) 

(1.2) 


where  are  known  constants.  The  problem  may  be  linearized  by  noting 

that,  in  the  infrared  region  and  at  terrestrial  temperatures,  the  variation 
of  B due  to  changes  in  v is  much  less  than  the  variation  due  to  T. 

The  variables  may  be  separated  by  first  computing  the  brightness 
temperature  T„(v)  from  (1.2): 

D 


Tg(v)  = B ^[V,N(V)] 
C^V 


Jin 


C^V-" 

N(V) 


+ 1 


(1.3) 


The  variables  may  then  be  separated  by  normalizing  the  measured  spectral 
radiance  to  the  black-body  equivalent  at  a reference  frequency  : 


N^(V)  ^ B[v^,Tg(v)] 


(1.4) 


Then  a normalized  version  of  (1.1)  may  be  solved: 


N^(V)  = B[v^,T(p^)]t(V,p^) 


x(p  ) 


Blv,,T(p)l  d.(p) 
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(1.5) 


In  this  form,  the  equation  is  linear  in 
has  been  found,  the  temperature  profile 
equation: 

T(p)  = b"^[v^,B^.(p)] 


Bj.(p)  = Bfv^,T(p)]. 
estimate  is  found  by 


Once  B (p) 
r 

using  Planck's 


(1.6) 


1.2.2  Minimum  Information  Method 

In  the  minimum  information  method,  (1.5)  is 
perturbed  form  by  using  the  initial  guessed  profile, 
guess  profile  by  T^(p),  we  have  from  (1.5): 


N^(V)  - B(v^,T^(p^): 


x(p  ) 

f 


B[V..T^(p)l^^(^-y-  dx(p) 


r o 


written  in  a 
Denoting  the  first 


(1.7) 


Then  it  is  assumed  that  T (p  ) = T (p  ) by  virtue  of  the  fact  that  accurate 

Os  s ■ 

surface  temperature  estimates  are  available  from  the  "window  channel"  mea- 
surements. With  this  assumption,  the  perturbed  equation  is 


AN(v)  = N (v)  - N (V) 
r o 


(1.8) 


where 

ABlv^,T(p)l  = B(u^,T(p)l  - B[v^,T^(p)]  (1.9) 

Now  by  using  numerical  quadratures,  (1.8)  can  be  written  in  the  form 


n . 
1 


A 


N'(V.  ) 


where 


= y r . . b.  ; i = 1,2 M 

11  1 
J 


(1.10) 


b.  = AB|v  ,T(p.)l 
J r J 


(1.11) 


and 

In 


r..  is  computed  as  a 
ij 

practice  x is  chosen 


function  of  x(p)  and 
to  be  proportional  to 


the  integration  rule. 

2/7  ...  . 

p , which  gives  essentially 


equally-spaced  intervals.  Then,  using  trapezoidal  integration: 


I^ij  ^i,j+l  ; j=l 

■^i.j-l  - ^i,j+l  ; l<j<L 


(1.12) 


where  T..  = P-  1®  ordered  monotonically  increasing  with  j. 

Thus  j=l  corresponds  to  the  top  of  the  atmosphere  (10  mb  in  this  study), 
while  j=L  denotes  the  1000  mb  pressure  level.  For  DMSP,  L = 70. 

It  is  convenient  to  write  (1.10)  in  vector  form  as 


n = Fb  (1.13) 

where  n is  an  M-dimenslonal  column  vector  of  radiance  variations,  b is  an 
l.-di.-nensional  column  vector  of  brightness  temperature  variations  and  F is  the 
matrix  of  transmittance  functions.  Since,  in  practice,  L > M,  the  problem  of 
inverting  (1.13)  to  find  b is  ill-conditioned  and  the  solution  is  non-unique. 
However,  once  a solution  for  b is  found,  the  temperature  profile  can  be  found 
by  using  (1.1),  (1.9),  and  (1.6). 

In  order  to  find  an  appropriate  value  of  b,  consider  the 
following  generalized  least  squares  problem.  We  wish  to  fi.-.d  b to 
minimize 

J(b)  = b^  E~^b  + (n-Fb)^  R~^(n-Fb)  (1.14) 


where  E and  R are  positive  definite  symmetric  matrices.  The  value  of 
b which  minimizes  J(b)  is 

b = ef'^CFef'^  + R)"^n  (1.15) 


This  equation  may  be  interpreted  in  terms  of  Bayesian  estimation  theory 
if  several  assumptions  are  made.  Assume  that  b is  a zero  mean  Gaussian 
random  variable  with  covariance  matrix  E.  Similarly  assume  that  n - Fb 
is  a zero-mean  Gaussian  random  variable  with  covariance  matrix  R.  Then 
b is  the  conditional  mean  of  b given  n and  is  also  the  minimum-variance 
estimate  of  b.  Eq.  (1.15)  is  also  in  the  form  of  the  measurement  update 
equation  for  the  celebrated  Kalman  Filter  (Kalman^  ^“^  ^).  in  practice 
the  values  of  E,  R and  F will  not  be  known  precisely  and  the  random 
variables  will  not  in  general  be  Gaussian.  Nevertheless,  this  technique 
can  yield  quite  good  performance.  In  the  vertical  temperature  retrieval 
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literature,  (l.l)  forms  the  basis  of  the  "foil  statistical"  method 
(cf.  Fritz  et  al.^  ^ and  is  generally  used  several  times  in  an  attempt 
to  iteratively  converge  to  the  correct  solution. 

The  minimum  information  method  uses  a special  case  of  (1.15) 


and  is  found  by  setting 

9 9 

E = a;:i.  R = o;i. 
h K 


(1.16) 


T T - 1 

b . = r (rr  + yn  n 

mi 


(1.17) 
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where  y = 0_/o'  is  the  "noise-to-s ignal " power  ratio.  This  solution  was 
K h 

considered  by  Foster  (1.15)  who  noted  that  it  corresponded  to  a maximum 
entropy  ensemble  in  which  there  is  a minimum  of  information  regarding  the 
statistical  characteristics  of  the  random  variables  b and  n.  This  is 
equivalent  to  the  assumption  of  isotropic  scattering  about  the  mean  for 
each  random  variable  and  further  that  the  random  variables,  under  the 
Gaussian  assumption,  are  independent. 

Twomey^  ^ ^investigated  problems  of  this  type  and  suggested 
the  use  of  a trial  solution.  The  final  solution  would  then  be  constructed 
to  minimize  the  deviation  from  the  trial  solution.  The  elements  of  this 
process  are  present  in  both  the  "full  statistical"  and  the  minimum-informa- 
tion solutions.  As  E 0 (or  y > m)  for  fixed  R,  b -»■  0 and  the  estimated 
temperature  profile  approaches  the  initial  guessed  profile.  This  also 
holds  if  R with  E held  fixed.  In  this  case  either  the  confidence  in 
the  initial  guessed  profile  is  high  or  the  measurements  are  extremely  bad. 

The  opposite  situation  occurs  when  confidence  in  the  initial  guess 
is  low  (E  -*■  CO,  with  R fixed)  or  the  measurements  are  extremely  good  (R  -*  0, 
with  E fixed).  In  either  case  Y 0 and  b is  obtained  via  the  pseudoinverse: 


^ T T -1 
b = EF  (FEP  ) n 


(1.18) 


Note  that  n = Fb,  thus  satisfying  (1.13),  and  in  addition  the 
weighted  cost 


T -1 

J(b)  = b E b 


(1.19) 


is  minimized. 


This  method  is  also  closely  related  to  the  method  of  Backus  and 

( 1 7 ) 

Gilbert  who  investigated  the  non-uniqueness  problems  associated  with 


r 


cietormiuat  ion  of  fine  geophysicai  striu-tiiro  from  gross  earth  data  measurements. 
In  their  metliod,  tl)ey  find  a solution  b which  minimizes  the  quadratic  cost 

(b-b^)*^  W~^(b-b^)  (1.20) 

where  b^  is  the  guessed  brightness  profile  and  W is  a positive-definite 
symmetric  weigliting  matrix.  They  also  constrain  the  estimated  measurements 
to  correspond  to  the  actual  measurements: 

h(b)  = n (1.21) 

% 

where  n denotes  the  vector  of  actual  radiance  measurements.  The  solution 
is  found  by  using  the  following  iteration; 

;(k-bi)  ^;(k) 

In  - h(h^^^)  - (b^  - a. 22) 

whore 

I 

j (1.23) 

3b  I 

iV"' 

Note  that  G^*^^  ^ I’  from  (1.13).  This  iteration  is  thus  similar  in  form 

to  the  "full  statistical"  approach  except  that  an  additional  term 

G^*^hb^  - b^'^h  is  present.  However,  the  fundamental  difference  is  that 

n(b)  n,  wliile  in  the  full  statistical  or  minimum- informat  ion  approach 

(1.21)  is  never  satisfied.  It  is  interesting  to  note  that  in  their  numerical 

^ (k) 

studies,  Backus  and  Gilbert  simplified  (1.22)  by  assuming  b^  = b . This 
is  precisely  what  is  done  in  the  iterated  "full  statistical"  and  minimum- 
information  methods,  namely,  at  each  step,  the  best  estimate  of  the  initial 
profile  is  set  to  the  estimate  on  the  previous  iteration.  With  this  assump- 
tion, the  method  of  Backus  and  Gilbert  reduces  to  (1.18)  under  the  assump- 
(k ) 

tions  W = E and  G = T.  It  is  Interesting  to  note  that  their  method  does 


not  explicitly  account  for  statistical  uncertainty  in  the  measurements  as 
do  the  "full  statistical"  and  minimum-information  methods. 


The  minimum-information  solution  is  seen  to  offer  a compromise 
between  using  the  initial  guess  profile  and  matching  the  guessed  and  measured 
radiances.  The  virtue  of  using  the  first  guess  profile  is  that  erroneous 
atniospher i c waves  which  may  be  generated  as  a result  of  ill-conditioning 

are  suppressed.  However,  this  smoothing  effect  can  lead  to  biases  in  the 

( 2 ) 

troposphere  under  partially  cloudy  conditions. 

The  DMSP  temperature  retrievals  are  performed  in  the  SOUNDER 
program,  in  wliich  up  to  five  iterations  of  (1.17)  are  made.  Using  the 
notation 

c = r^d’r'^  + Yi)~^ 

the  iterations  proceed  as  follows.  On  the  k*'^'  iteration 


„"(k) 


1 + C.  n - n I 


where  T is  the  brightness  temperature  computed  from  the  temperature 

profile  via  Planck's  equation,  h is  the  measured  spectral  radiance  vector, 

^(k-1)  *(k-l) 

and  n is  the  predicted  radiance  found  from  T and  the  radiative 

transfer  model.  Tlie  surface  temperature  estimate  is  not  iterated.  The 

iteration  terminates  wlien  the  radiance  deviation  in  each  channel  falls 

below  a preset  thresiiold.  In  DMSP,  all  thresholds  are  set  equal.  It  should 

be  pointed  out  that  tliis  may  lead  to  problems,  since  the  uncertainty 

associated  with  each  channel  is  not  the  same.  Inspection  of  individual 

iterations  reveals  that  channels  1,2,3  are  almost  invariably  the  last  to 

meet  the  threshold  requirement,  which  is  in  intuitive  agreement  with  the  fact 

that  expected  errors  for  these  channels  are  greater.  This  may  induce  errors 

larger  than  necessary  in  the  temperature  estimates  above  200  mb. 


1.2.3  Regression  Analysis 

The  use  of  regression  analysis  is  a powerful  tool  in  analyzing 
the  relationships  between  sets  of  variables  of  different  types.  In  the 
present  problem,  we  are  interested  in  analyzing  the  relationship  of  vertical 
temperature  profiles,  actual  measured  spectral  radiances,  initial  forecast 
fields  (first  guess  fields)  id  the  TTPACK  analysis. 

The  several  objectives  of  the  analysis  were: 

(1)  Investigate  the  information  content  of  tlie  initial 
foreiast  field  and  spectral  radiances  relative  to 
vertical  temperature  profiles. 


(2)  Perform  a validation  study  of  the  minimum-information 
technique  relative  to  an  optimum  linear  processor. 

(3)  Develop  methods  for  improving  the  minimum- in format ion 
technique . 

(4)  Perform  analyses  to  determine  robustness  of  estimators. 


It  was  decided  that  regression  analysis  would  be  used  for  developing 
statistical  models  rather  than  principal  components  (statistical  eigenvector 
methods)  for  several  reasons.  First,  a weakness  of  the  principal  components 
methods  is  that  it  depends  on  the  units  of  measurement.  If  a variable  is 
measured  in  such  small  units  that  its  numerical  values  dominate  those  of  the 
other  variables,  the  first  principal  component  will  essentially  reflect  the 
behavior  of  this  variable.  This  is  sometimes  handled  by  normalizing  each 
variable  to  zero  mean  and  unit  variance.  However,  the  correct  weighting  to 
be  given  to  each  variable  is  not  clear,  since  it  depends  on  the  nature  of 
the  problem.  The  most  significant  variables  might  be  given  higher 
weighting.  In  the  sounding  problem,  this  is  equivalent  to  saying  that 
accurate  temperature  retrievals  may  be  more  important  in  certain  pressure 
regions  than  in  others.  Secondly,  by  working  in  the  natural  units  of  the 
problem,  the  results  are  made  easier  to  interpret.  Thirdly,  the  first  few 
principal  components  computed  by  Smith  and  Woolf ^ ^ ^ and  Smith ^ ^looked 
sufficiently  similar  to  the  transmission  functions  to  indicate  that  regression 
models  may  be  competitive  in  practice.  Finally,  the  objectives  of  this 
study  were  limited  to  investigation  of  the  minimum-information  solution  and 
methods  of  improvement  rather  than  investigation  of  tlie  many  possible 
alternatives  to  temperature  retrievals. 

1 . 2 . 3 . 1 Regression  Models 

Two  different  regression  models  were  utilized  in  an  attempt  to 
assess  the  performance  of  the  minimum  information  solution.  For  temperature, 
the  estimate  for  the  i—  level  is 
It  8 

i=l  ^ k=l 


T.  = ) a.  .T  . + )■  B.,  r,  + Y,  ; i 

1 ij  0l  , ik  k 'j 


1,....13 


.th 


where  T^^  is  the  first  guess  at  the  j--  level  and  rj^  is  the  measured  spectral 


th 


radiance  for  the  k—  ciiannel.  The  coefficients  a..,  B.,  and  Y.  are  determined 

ij  ik  i 

by  means  of  a stepwise  regression  tecimique  in  which  only  the  most  significant 


coefficients  are  retained  in  the  model.  The  units  for  temperature  are  °C 

2 


and  the  units  for  the  radiances  are  ergs  (cm 

For  estimating  lieights,  the  model  is 

n 8 


-1,-1 

sec  - ster  - cm  ) 


H. 
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y 6 . . H . 

,i=i 


L ^ -1,  T b • • • -i  ' 

ik  k 1 


k=l 


.th 


where  . is  the  first  guess  height  (in  ft)  at  the  J—  level.  The  regression 


coefficients  <5, 


and  are  found  by  tlie  same  stepwise  regression 


method  used  for  the  temperature  estimates.  The  regression  method  used 
was  performed  by  the  program  RLSEP,  which  was  available  at  AFGL  as  part 
of  the  International  Mathematical  Statistics  Library  (IMSL) . 

It  is  important  to  point  out  that  the  regression  method  implicitly 

takes  into  account  empirical  linear  correction  techniques,  such  as  the  method 

( 7 ) 

of  Weinreb  and  Fleming.  In  this  method  the  observed  systematic 

disagreeement  between  temperature  inversions  and  RAOBs  is  reduced  by  employing 
a linear  correction  to  the  spectral  radiances.  With  the  usual  relation 
between  brightness  temperature  and  radiances  given  as 
b = Cn, 

Weinreb  and  Fleming  proposed  adding  a term  to  n to  give 
b'  = C(n  + 6n) 


The  value  of  bn  was  then  found  using  a least  squares  fit  to  a set  of 
actu.al  RAOB  data.  This  method  is  used  in  the  AFGWC  SOUNDER  program.  The  necessity 
of  using  such  a technique  is  an  indication  of  the  use  of  statistical  methods 
in  present  use  and  of  the  usefulness  of  statistical  techniques  in  general. 

1 . 2 . 3 . 2 Stepwise  Regression 

In  computing  linear  regression  models,  it  is  important  to  deter- 
mine the  effect  of  using  different  subsets  of  independent  variables.  It  is 
desirable  to  use  only  those  variables  which  significantly  reduce  the  regres- 
sion errors.  Overfitting  can  lead  to  models  which  are  overly  sensitive  to 
noise.  Stepwise  regression  provides  a partial  automation  of  the  variable 
selection  process.  It  is  based  on  a technique  which  in  the  process  of 
computing  an  ordinary  regression  on  m predictors  obtains,  at  essentially  no 
extra  expense,  m intermediate  regressions  which  are  useful  in  determining 
functional  relationships  between  the  dependent  variable  (temperature  or  height) 


and  several  selected  subsets  of  the  total  set  of  independent  variables, 
or  predictors.  In  the  simplest  case,  one  variable  is  added  at  each  step. 

There  are  several  meaningful  statistical  criteria  which  may  be  used  to 
determine  which  predictor  variable  to  add  to  the  model: 

(i)  add  the  predictor  whose  F-to-enter  statistic  has  the 
largest  value. 

(ii)  add  the  predictor  that  gives  the  greatest  decrease  in  the 
residual  sum  of  squares. 

(iii)  add  the  predictor  which  gives  the  greatest  increase  in  the 
multiple  correlation  between  the  dependent  variable  and  the 
predictors. 

The  F-to-enter  statistic  for  a predictor  is  the  F-statistic  for 

testing  the  significance  of  the  regression  coefficient  the  predictor  would 

have  if  it  were  added.  All  four  of  these  criteria  are,  in  fact 

mathematically  equivalent;  however  the  F-to-enter  statistic  was  used  in 

programs  utili.'^ed  in  this  study. 

In  addition  to  adding  variables  it  is  also  desirable  to  provide 

for  removal  of  variables.  It  may  happen  that  recently  added  variables 

may  in  combination  render  a variable  entered  at  an  earlier  state  statistically 

insignificant.  Removal  of  variables  is  based  on  an  F-to-remove  test. 

In  the  stepwise  regression  procedure,  the  test  for  removal  of  predictors  is 

first  made  for  each  predictor  already  in  the  model.  If  a predictor  is  removed, 

the  program  proceeds  to  the  next  step.  If  no  predictor  is  removed,  the  step 

continues  with  the  addition  of  a predictor.  The  stepwise  regression 

procedure  terminates  when  no  predictor  is  either  deleted  or  added  at  a 

step.  The  details  of  the  stepwise  regression  procedure  may  be  found,  for 

(18) 

example,  in  Draper  and  Smith.  ° 

1,2.4  Discussion  of  DMSP  Analysis 

A simplified  information  flow  diagram  for  the  DMSP  minimum- 
information  technique  is  shown  in  Fig.  1.1,  depicting  the  major  inputs, 
outputs,  computing  blocks  and  parameters. 

As  a result  of  investigating  the  minimum-information  techniques 
and  discussions  with  AFRI.  and  AFGWC  personnel,  several  possible  sources  of 
error  have  been  defined.  In  approximate  order  of  decreasing  significance 
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they  are ; 


(1)  Clouds . Currently,  the  DMSP  software  selects  a 
single  candidate  sounding  from  each  240  nautical  mile 
grid  point  to  generate  the  clear  column  radiance  file. 

This  candidate  is  chosen  purely  on  the  basis  of  the  ratio 
of  the  observed  radiance  of  an  upper  air  channel  to  the 
window  channel.  While  this  test  is  employed  to  eliminate 
cloudy  areas,  discussions  with  AFGL  personnel  indicated  that 
considerable  cloud  contamination  may  occur.  Future  work  will 
be  necessary  to  define  and  evaluate  cloud  contaminant  reduc- 
tion techniques,  such  as  the  adjacent  field  method. 

(2)  Water  Vapor  Correction.  The  transmittance  functions  of 
the  atmosphere  are  not  updated  for  the  contribution  of 
water  vapor  arising  from  changes  in  the  temperature  of  the 
profile  (via  a constant  mixing  ratio  assumption)  as  the  algorithm 
proceeds . 

(3)  CO-  Transmittance.  Errors  in  cause  errors  in  the 

I CO2 

inversion  process  since  this  parameter  is  used  directly 
in  the  inversion.  At  present,  no  temperature  dependence  is 
assumed . 

(4)  Temperature  Profile  Adjustment.  The  principal  sources  of 
error  here  appear  to  be  the  problem  of  determining  the  trans- 
mittance matrix  F accurately  and  defining  a more  complete 
error  model. 

(5)  Empirical  Corrections.  This  technique  is  used  to  correct 
for  bias  before  the  inversion  is  performed. 

(6)  Channel  Width/Frequency.  Errors  in  assumed  channel  center 
frequency  and  width  cause  errors  in  the  computed  temperature 
profile. 

Other  possible  sources  of  error  include  level  conversion, 
layer  spreading,  and  earth  location. 

The  minimum-information  solution  which  is  being  used  for  temperature 
inversion  is  a prime  candidate  for  a source  of  a significant  percentage  of  the 
observed  errors.  One  of  the  assumptions  employed  in  the  model  is  that  the 


31 


errors  associated  with  each  IR  channel  have  equal  variances,  (cf.  Section 
1.2.2).  The  accuracy  of  this  assumption  is  questionable.  Another  assump- 
tion, which  is  more  questionable,  is  that  the  temperature  estimation  errors 
associated  with  each  pressure  level  are  of  equal  variance.  It  is  further 
assumed  in  the  model  that  all  errors  (both  radiance  measurement  and  tempera- 
ture inversion)  are  uncorrelated.  This  assumption  does  not  appear  to  be  met 
in  practice. 

The  transmittance  matrix  F is  formed  by  taking  differences  of  the 
transmission  functions,  so  that  errors  in  the  transmission  functions  translate 
directly  in  errors  in  the  transmittance  matrix.  One  source  of  error  is  the 
assumption  that  F is  independent  of  temperature.  Also,  F is  determined  using 
models  of  water  vapor  transmittance,  ozone  transmittance  and  CO^  transmittance 
and  thus  is  dependent  on  accuracy  of  all  of  the  models. 

The  minimum-information  solution  may  be  modified  by  doing  several 
things.  First,  the  adequacy  of  the  noise  models  should  be  validated.  This 
has  been  done  by  performing  statistical  analyses  on  the  data  available  at 
AFGL.  Covariance  matrices  for  both  the  measurements  and  estimation  errors 
have  been  estimated. 

The  matrix  F can  also  be  adjusted  to  achieve  better  performance. 

Two  criteria  can  be  used;  (1)  variance  matching,  (2)  whiteness  tests. 

Variance  matching  is  easier  to  perform,  but  whitening  the  residuals  (differ- 
ence between  measurement  (radiance)  and  estimated  measurement  (brightness 
temperature),  is  a more  powerful  approach.  A hybrid  approach  which  is 
similar  to  whitening  the  residuals  is  used  in  the  sequel  to  study  the  proper- 
ties of  the  minimum-information  method. 

Another  problem  with  the  minimum-information  method  is  that  it 
is,  in  its  present  form,  an  ad  hoc  approach.  From  a Bayesian  point  of 
view,  only  a single  iteration  need  be  taken  since  the  model  is  linear  with 
no  dynamics  and  a single  vector  measurement.  This  problem  has  been  addressed 
implicitly  in  this  study  by  comparing  minimum-information  retrievals  with 
statistical  non-iterative  methods. 
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2 MINIMUM  INFORMATION  TECHNIQUE  AT  AFGWC 


The  minimum  information  technique  is  part  of  an  AFGWC  soft- 
ware system  called  E- package  or  EPKG.  Most  of  the  subprogram 
components  of  the  system  are  written  in  the  language  of  FORTRAN  V 
and  in  the  particular  dialect  which  is  resident  on  the  UNIVAC  1100  Series 
Operating  System.  Extensive  use  has  been  made  of  the  unique  features 
of  FORTRAN  V and  1100  Operating  System  in  the  EPKG.  Because  of  this, 
implementation  of  the  EPKG  on  other  systems  would  be  a very  involved 
process  and  difficult,  particularly  if  undertaken  by  one  not  familiar  with 
the  UNIVAC  Operating  System  and  FORTRAN  V. 

The  procedures,  calculations,  etc.  which  constitute  the  mini- 
mum information  technique  are  embodied  in  the  main  program  SOUNDER 
and  its  associated  subprograms  of  the  EPKG.  An  annotated  flow  diagram 
showing  the  major  procedural  and  computational  steps  in  SOUNDER  is 
given  in  section  2.1.  The  functions  of  the  simpler  subprograms  called 
by  SOUNDER  are  described  on  the  flow  diagram.  The  functions  of  the 
more  complex  subprograms  are  described  in  sub-sections  of  2.1  verbally, 
diagrammatically  or  by  a combination  of  these  forms. 


33 


2.1  MAIN  PROGRAM  SOUNDER 


The  following  annotated  flow  diagram  shows  the  major  pro- 
cedural and  computational  steps  in  SOUNDER  which  relate  to  the  im- 
plementation of  the  minimum  information  technique.  The  operations  of 
SOUNDER  which  deal  with  the  locating  or  addressing  and  storage  and 
retrieval  of  data  from  the  AFGWC  data  bases  are  not  shown.  The 
decision  criteria  described  as  well  as  the  values  quoted  for  constants 
and  test  variables  were  those  in  use  in  the  operational  version  at  AFGWC 
as  of  1 April  1977. 

Various  levels  of  detail  have  been  used  in  describing  the 
operations  in  SOUNDER  and  its  subprograms.  The  general  guidance  used 
in  determining  the  level  of  description  detail  was  the  specificity  and/or 
uniqueness  of  the  operation  to  the  implementation  of  the  minimum  in- 
formation technique  at  AFGWC.  For  example,  subprogram  elements 
ITERAT,  H20TAU,  and  COFWGT  are  described  in  considerable  detail  as 
they  constitute  the  essence  of  the  minimum  information  technique.  On 
the  other  hand,  subprogram  elements  THKNES,  MANTMF,  and  TRPFIN 
which  employ  techniques  familiar  to  most  meteorologists  for  the  deter- 
mination of  thicknesses,  the  height  and  temperatures  of  mandatory 
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levels  and  the  location  of  the  tropopause  from  lapse  rate  changes  are 
described  very  briefly.  In  between  these  two  extremes  are  subprogram 
elements  INITMP  and  SMOOTH  which,  although  for  the  most  part  use 
conventional  meteorological  procedures,  also  have  a few  somewhat 
unique  but  critical  functions  in  the  implementation  of  the  minimum  in- 
formation technique. 
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Subroutine  THROW  performs  the  following 
unctions: 

1 - Terminates  analysis  of  a sounding 

whose  underlying  surface  was  land 

2 - Terminates  analysis  of  a sounding 

whose  zenith  angle  exceeds  60° 

3 - Determines  zenith  angle  index  number 

of  stored  CO^  transmission  function 


Subroutine  TRAMIS  determines  the  drum 
sector  addresses  of  the  CO2  and  O3 
transmission  functions  for  the  satellite, 
month  of  the  year,  latitude  and  zenith 
angle  of  the  sounding  and  then  inputs  these 
fuuctvons  from  mass  storage.  A description 
of  the  stored  transmission  functions  is 
given  in  Table  2 . 2 


Subroutine  INITMP  uses  stored  climato- 
logical profiles  and  AFX3WC  analysis 
temperature  fields  to  construct  first  guess 
profiles.  Subroutine  SMOOTH  smoothens 
the  transition  from  the  climatologically 
derived  and  analysis  sections  of  the 
first  guess  profile. 
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Suhroutmo  HERAT  executes  ihc  specific 
procedures  of  the  nunimuni  information 
technique  to  retrieve  the  temperature 
profiles . 

Subroutine  H20TAU  calculates  the  com- 
bined carbon  dioxide-,  water  vaix')r  and 
o^ont  transmission  function  for  channels 
1 through  7. 


Subroutine  COPVVGT  generates  the 

coefficient  matrix 

sec  Ref.  2 equation  12. 


Subr.’ut.ru;  TliKNFS  • 
in ic< nesses  eetwfe:. 
tr.-  •.inaat.iy  bvel.. 

I etrii  vovi  *em;  t r ilure 
assumed  water  vaiv  r 


liculates  tin 
i:ni  hei.jnts 
ttn 

pr.dile  ana 
profile. 


‘•t 


an 


maktmp 


J 


t 


I THPriN 


r . 

Stop  ^ 

/ 

- . ^ 


Subroutine  CMTEMP  retrieves  the 
temperature  profile  by  the 
statistical  C matrix  solution 
technique. 


Subroutine  MANTMP  calculates 
temperatures  at  the  mandatory 
levels  by  interpolation  between 
the  computational  levels. 


:v) 


Subroutine  TRPPIN  determines  the 
tro^.'Opause  level  in  the  final 
temperature  profile. 


Table  2.1  Data  Items  in  a Radiance  Set  of  the  READYRADIANC  File 


Name 

Description 

Latitude  Belt  No. 

An  integer  from  1 to  45  denoting  the 

4 degree  latitude  belt  in  which  the 
sounding  was  made.  North  Pole  is  belt 
no.  1,  South  Pole  is  belt  no.  45 

Longitude 

Longitude  of  sounding 

Latitude 

Latitude  of  sounding 

Minutes 

Minutes  after  the  hour  of  sounding 

1 

Hemispheric  Indicator 

Indicator  which  shows  the  hemisphere  of 
the  sounding 

3D-NEPH  Location 

The  3D-NEPH  box  no,  I and  J coordinates 
of  sounding 

Julian  Hour 

Julian  hour  of  the  sounding 

Channel  Radiances 

Measured  radiances  for  the  eight  channels 

Zenith  Anglo 

Zenith  angle  for  the  measured  radiances 

Terrain  Height 

Terrain  height  at  the  location  of  the 
sounding 

Surface  Indicator 

Indicator  denoting  the  type  of  surface 
(sea,  land,  ice)  underlying  the  sounding 

Mandatory  Level 

Analysis  Data 

AFGWe  analysis  temperature  and  D values 
for  the  fifteen  mandatory  levels  from 

1 000  mb . to  1 0 mb. 
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Name 


Description 


Tropopause  Data 

AFGWC  analysis  tropopause  temperature 

and  pressure 

Dummy  Block  Station 

A number  which  identifies  the  satellite 

Number 

from  which  the  sounding  was  made  and 

additionally  denotes  information  relevant 

to  the  filing  of  the  sounding  data 

r 
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Tabic  2,2  Description  of  the  Stored  Carbon  Dioxide  and  Ozone 
Transmission  Functions 


Carbon  Dioxide 

Ozone 

No.  of  atmospheric  pressure  levels  at 
which  transmission  from  the  satellite 

IS  specified 

70 

70 

No.  of  satellites 

4 

1 * 

No.  of  channels 

6 

6 

No.  of  atmospheres 

6 

5 

No.  of  zenith  angles 

15 

1 ** 

* Ozone  transmissions  are  not  satellite  specific 
**  Ozone  transmission  is  specified  for  0°  zenith  angle  only 


Satellites  - WX6530,  WX8531,  WX9532,  WX1534 

Carbon  Dioxide  Atmospheres  - 15°  N Annual,  30°  N Winter 

30°  N Summer,  U.S.  Standard 
60°  N Winter,  60°  N Summer 

Ozone  Atmospheres  - Tropical,  Mid-latitude  summer. 

Mid-latitude  winter.  High  latitude 
summer.  High  latitude  winter 


J . 1 . 1 Subprogram  Klrnicnt  INITMP 


Subroutino  INITMP  builds  the  first  guess  temperature  profile. 

The  temperatures  and  heights  are  specified  at  70  computational  levels 
whicli  are  related  to  pressure  by  the  equation 

N ■ apV7  + ^ 1) 

wliere 

N = computational  level  number  (integer) 

P - pressure  (mb) 
a = 9.9587402862 
8^-1  .b7  16270699 

The  constants  a and  8 ore  determined  by  requiring  that 
N = 1 when  P = .01  mb. 
and  N = 70  when  P = 1000  mb. 

For  the  computational  levels  above  10  mb.  , temperatures  corres- 
ponding to  one  of  six  climatological  atmospheres  listed  in  Table  2.3 
arc  used.  The  particular  climatological  profile  selected  is  determined 
by  the  month  and  latitude  of  the  sounding. 

For  the  computational  levels  which  correspond  to  the  mandatory 
levels  from  10  mb.  down  to  1000  mb.,  the  ARTWC  analysis  tempera- 
ture's and  D values  are  used  for  the  first  guess  temperatures  and  heights. 
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Tabic  2.3  Climatological  Temperature  Profiles  Used  Above  10  mb. 


1 . 

15°  N 

or  S Annual 

2. 

30°  N 

Winter 

3. 

30°  N 

Summer 

-1. 

45°  N 

Spring,  Fall 

5. 

60°  N 

Winter 

b. 

o 

O 

Summer 

If  a superadiabatic  lapse  rate  results  in  850  to  1000  mb.  layer, 
the  temperature  at  1000  mb.  is  reduced  by  1°C  while  the  temperature 
at  850  mb.  is  increased  by  .5°C.  Successive  replications  of  this 
procedure  are  made  until  the  lapse  rate  in  the  layer  is  no  longer 
superadiabatic . 

INITMP  calculates  the  temperatures  at  the  computational  levels 
between  the  mandatory  levels.  The  interpolation  algorithm  for  the 
temperature  at  a computational  level  lying  between  any  two  mandatory 
levels  which  do  not  bracket  the  tropopausc  is  given  by 
Ti  = T,^  - (Tn,  - Tni-i)  * Ri 

where  _ In  (Pj)  - In  (Pm-i)  (2.2) 

In  (Ppi)  - In  (Ppi-i) 
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Tj,  Pi  = temperature,  pressure  at  the  computational 
level 

^m < ^m  ~ temperature,  pressure  at  the  first  mandatory 
level  below  the  computational  level 

'^m-1'  ^m-l  ~ temperature,  pressure  at  the  first 

mandatory  level  above  the  computa- 
tional level 


The  interpolation  algorithms  used  when  the  mandatory  levels  m and  m-1 


bracket  the  tropopause  are  given  by 

Ti  = Tm_i  + (Tt  - Tm-l)  * ~ (Pm- 1) 

In  (Pt)  - In  (Pm-l) 


(2.3) 


when  the  computational  levels  lie  above  the  tropopause  and 


In  iP^)  - In  (Pj 

^i  ^ - T^)  * 


(2.4) 


when  the  computational  levels  lie  below  the  tropopause.  Tt  and  P^  are 
the  temperature  and  pressure  of  the  tropopause. 


2.1.2  Subprogram  Clement  SMOOTH 

Subroutine  SMOOTH  passes  a three  point  smoother  thru  the  first 
guess  temperature  profile  from  level  21  upward  to  level  10.  The 
purpose  of  this  is  to  remove  any  discontinuity  that  might  exist  between 
that  portion  of  the  first  guess  temperature  profile  obtained  by  interpo- 
lation of  the  AFGWC  mandatory  level  temperatures  and  the  climatological 
first  guess  data  used  above  10  mb.  The  smoothing  algorithm  is  given 
by 

Ti  - * T,_i  ^ ^2.,  * T,  + W3,i  * Tj  + i (2.5) 

The  weights  Wj  ^ are  given  in  Table  2.4. 


Table  2.4  First  Guess  Temperature  Profile  Smoother  Weights 


2.1.3 


Subprogram  Element  HERAT 


CtilculaK'  tr.insnubSioii  functions  dut-  to 
wjttT  v.ipor,  ozone  and  carbon  ihoxide 
at  computational  levels  I throuqh  70  for 
channels  1 through  7. 


If  abnormal  rc-turn  tror..  1120TAI\  execute 
abnorn:al  return  to  SOUNDI'R. 


Compute  tne  cot  fficier.l  maliix^ 
{Equation  12,  Ref.  2 ) 


Planck  lunction  for  the  surface  level 
and  computational  level  70  (iOOn  mb 
level)  arc  equal. 


Calculate  equivalent  radiative  temperature 
tor  each  measured  radiance,  then  calculate’ 
Planck  radnneo  tor  those  tcmpeiaturos  it 
tne  reference  frequency. 
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Where 

U>K(1)  = 

DRBARdl 
SR/\n(J)  = 
RBAHd)  = 


C‘)rrection  to  be  added  to 
I’qujvjlont  measured  radiance 
of  ehann-'l  1 
” D ' matrix 

Mean  radiance  correction 
vector 

Pv^uivalent  measuied  raaiancr 

0}  channel  I 

Mean  radiance  vectiii 


DFLViT^D  ^ SRAD{KR)-  BSrC(NR)-CR^\D(NR) 
where 

SRAD(NR)  = "D  " n-.dtrix  corrected 

equivalent  measured  r.tdiqncc 
of  window  channel  NR 
BSFC(NR)  ,CRAD(NR)described  pret’iously 


It  the  calculated  window  ctiannel  diffi'rence 
IS  less  than  -CLDDlf,  the  sounding  is 
assumed  to  he  too  cloud  contaminated  near 
tne  suttace  lot  tuithei  analysis  and  an 
abnormal  leturn  to  SOUNDFR  is  executed. 
CLDDIF'  is  currently  set  equal  to  1.1. 


■ HI 


r 


I 


If  calculated  window  channel  difference  is  greater 
than  HOTSrC,  execute  an  abnormal  return  to 
SOUNDER  as  the  sounding  cannot  be  processed 
further  either  because  the  underlying  surface 
temperature  is  too  hot  or  the  initial  temperature 
profile  guess  is  bad.  HOTSFC  is  currently  sc-t 
equal  to  1 . 1 


It  the  sum  squares  of  tne  deviations  SMSQ  is 
less  than  or  equal  to  SNOIS  ana  the  maximum 
deviation  DIW.AX  is  less  than  EPS,  jump  to  15 
and  compute  final  temperature  profile,  SNOIS 
is  current'y  equal  to  1 . 5 and  EPS  equal  to  . 5 . 


Ai^normal  return  to  SOUNDER  if  mviximum  number 
of  Iterations  have  been  performed. 


The  coefficient  n.atrix  ^ computed  in  COFWGT  is 
used  on  the  vector  of  deviations  of  calculated 
from  measured  radiances  to  compute  a new  Planck 
Function  Profile. 
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2.1.4  Subprogram  Element  H20TAU 


Iniry  Poii-t  H20TAI: 


SuhrouliJit'  tSlIjiS 

s th»'  situr^tion  vajv)r  pr'^ssure, 
.ivcriga  lenipCTdturtJ,  and  associatc'd 
tt^ibperjture  functions  for  the  "water 
vapor  layers".  The  "water  vapor  layers" 
start  at  level  43  anvi  end  at  level  70, 


Using  a fixed  relative  humidity  profile  for 
the  water  vapor  layers  (See  Table  2.5) 
ealculalc  water  vapor  pressure  in  each 
layt  I and  the  cjrnount  of  .vaier  vapor 
traversed  in  the  layers  givt-.n  the  .tenitn 
mole  Of  the  sounding. 


Calculate  tor  i acn  .'f  the  water  va^^u 
Ij/i.rs  the  pnssurc,  temperature  and 
-.vatei  vipt^r  } ath  lenoth  dependent  terms 
o:  the  witn  vapor  ti  insmission  functions. 


Calculate  tru  water  v-nvir  transmission 
t'.v  the  window  channel,  i.o,  channel  7, 
at  C'^i  putaticnal  levels  44  inrough  70. 
Levels  1 through  4 3 set  transmission 
equal  to  . 


lnit.alii-(  umsmission  for  channels  4.  5 
ina  I at  ill  computational  levels  to  I . 
Cilculatt  water  vainar  transmission  at 
coe  pul  Ui  Mtal  levels  4B  through  70. 


Test  cak  elated  witer  wajxir  Iransnussions . 
If  ttansnnssion  at  any  level  is  negative, 
execute  abnormal  return  to  ITLPjM, 


t 

Initi  ili/.o 
Trciiisnubsicnb  loi 
Channels  l , ^ ^ } 


Initialize  transmissions  for  channels  1 , 
2 & 3 at  all  computational  levels  to 
CO2  transmissions . 


Initialize 
Transmissions  at 
Levels  1 through  -13 
for  Channels  -1,  5 <Sc  6 


j C ilculati  Intermediate 
, iiansmission  Values  at 
I Levels  44  inrouq:;  70 
i for  Cnannels  4,  5 ^ c 


Calculate  final 
Transn-.issior.  V>alucs  ■ 
Ail  Levels,  Channels' 
i through  u 


Initialize  iransn.issions  for  channels  4, 
5 & h at  con.puiational  levels  1 through 
43  to  CO2  transmission. 


Calculate  ir.iermeduite  transmission 
values  equal  to  product  of  CO2  and 
H2O  transmissions  at  levels  44 
through  70  tor  channels  4,  5 ^ n. 


Calculate  ri<AN , final  transmission 
I'i?AN  TiUNI  » TK>\NO  ■'’fCO 

TRAM  - initial  or  intermediate  value 
TR/\N(')  - ozone  transmission 
^ = zenith  angle 


It  trinsmission  at  any  level  tor  .iny 
Channel  is  negative,  execute  adnormal 
retuin  to  ITfKAT. 
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Tabic  2.5 

Relative  Humidity  Profile  in  tlie  Water  Vapor  Layers 


'computational  ; Average  Relative  ' Computational  ■ Average  ] Relative 
I Level  No.  j Layer  I Humidity  , Level  No.  i Layer  ! Humidity 
I ^ Pressure  : , Pressure  ' 


2.1.5 


Subprogram  Element  COTWC-T 
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2.2  GAC  RESEARCH  VERSION  OF  SOUNDER 

The  Geo-Atmosphcrics  Corporation's  (GAC)  Research  Version  of 
SOUNDER  is  a modification  of  the  SOUNDER  program  which  was  opera- 
tional at  AFGWC  on  1 April  1977.  Implementation  of  the  minimum 
information  technique  in  the  research  version  is  identical  to  that  in  the 
operational  version.  Major  modifications  that  were  made  involved  the 
input/output  of  data.  Exercising  GAC's  research  SOUNDER  provides 
information  attainable  on  weaknesses  and  strengths  of  the  operational 


version  of  SOUNDER.  j 

I 

The  GAC  research  version  like  the  operational  version  is  | 

written  in  FORTRAN  V and  must  be  run  on  a UNIVAC  1100  Senes  System. 

i 

The  program  requires  44K  words  decimal  to  load.  The  three  data  files 
required  are  the  carbon  dioxide  transmissions,  ozone  transmissions  and 
the  EPKG  statistics.  A new  main  program  was  written  which  generates 
a READYRADIANC  file  from  user  supplied  radiances,  first  guess  tempera- 
tures, and  heights. 

Additional  outputs  were  incorporated  into  SOUNDER  in  the  re- 
search version.  They  include  the  normalized  measured  radiances, 
radiances  calculated  from  the  guess  profiles,  D matrix  corrected  radiances, 
and  various  differences  in  these  quantities. 

1 

J 


Tho  sinqlo  qronU'st  value  in  tlie  C'.AC:  rt'search  version  of 
SOUNDER,  however,  lies  not  so  much  in  the  revised  proqram  itself, 
but  rather  in  the  research  capabilities  built  into  the  operation  and 
functions  of  various  elements  of  the  program.  Research  modifications 
have  enabled  one  to  easily  modify  any  of  the  preselected  data  or 
operations  incorporated  in  SOUNDEIf  such  as  the  assumed  water  vapor 
profiles  and  the  CLDDIE,  HOTSEC,  SNOlS  and  EPS  test  values. 

Other  more  complex  modifications  such  as  the  decoupling  of 
underlying  surface  temperatures  from  the  1000  mb  temperature  and 
provisions  for  forecast  data  inputs  allow  considerable  flexibility  in 
deriving  operational  sensitivity  and  performance  relative  to  available 
data  inputs. 

The  GAO  research  version  of  SOUNIhER  is  a powerful  tool  for 
exploring  hardware  and  software  deficiencies  and  improvements  upon 
routine  operational  performance. 


r 


3 ^NALYTR'.Al,  IIRKOU  ANALYSIS 


3.1  DATA  BASES 

3.1.1  Gent  r.il 

Three  data  btisos  used  in  this  study  were  derived  from  satellite 
infrared  radiometi'i  s , radiosoiKie  observations  (liAOb) , anti  Air  Foice 
Global  Weather  Centc'r's  (AEGWG)  13-hour  nunieriCMl  forec.ist  field.  In 
operational  practice,  the  toia’cast  data  are  useii  to  tabtain  a first-guess 
temperature  protile  foi  use  witii  tlu'  multi-channel  radiometric  uata  to 
obtain  a retrieved  lemperaturi’  piolile  asinu  tiie  minimum  intoimation 
technique.  This  proceouri'  is  oesci  iLiCtl  in  a iLdtiiled  sti  p by  stc'p  manner 
in  section  1 of  this  report.  lUVOB  datti  wen  colKcted  tor  only  those 
locations  wiu're  a I'Kar  (clouti-lia'e)  SLitelliti.-  souni.iinq  was  possible 
within  bO  nautical  miK'S  and  one  hour  in  time  of  one  anotiur.  All  satellite 
data  were  collected  lor  clouii-free  oci  an  reqions  only.  Data  trom  two 
satellite  birds  were  av<iilable,  Bird  9533  aiui  8531  , a morning  and  noon 
satellite , n'S pect i ve  1 y . 

Chiiracterist ics  associated  with  the  sati  llite  si  nsors  are  given 
in  Table  3.1.  Six  ol  tlu-  iiuiiometric  chaiiiu  Is  occur  on  the  |)eak  (15.  Ou  '”) 
or  wings  of  the  carbon  dioxide  molecule  abstirption- emission  band.  It  can 
be  seen  that  the  higlier  the  cliannei  number  tiie  greater  the  depth  ol 
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penetration  or  interrogation  of  atmospheric  conditions.  This  is  shown 
graphically  in  Fig.  3.1  where  the  distribution  of  the  weighting  function 
for  each  of  the  first  six  sensor  channels  shows  that  portion  of  the 
atmosphere  that  contributes  to  the  radiometric  power  received  at  the 
satellite. 

Channel  1 sensor  is  centered  on  the  peak  of  the  CO2  band  at 
15  n m or  668.5  cm~^  wave  number  and  receives  radiation  throughout  a 
half  power  band  width  of  3.5  cm  ^ with  a 0.25mw/m  • ster  • cm  ^ 
standard  deviation  of  the  nominal  noise  equivalent  spectral  radiances 
(NESN).  Channel  1 receives  upwelling  radiation  from  the  top  portions  of 
the  atmosphere,  centered  at  about  the  3i)  mb  level.  The  band  width  is 
increased  and  the  NESN  is  decreased  for  all  other  channels.  Channel  6 
weighting  function  peaks  at  the  1000  mb  pressure  level  whereas  channel 
7 (IZ.O^m  or  835.0  cm  M operates  in  a "window"  region  and  receives 
maximum  radiation  from  the  underlying  surface.  Channel  R operates  on 
the  18.7^  m or  535.0  cm”^  water  vapor  line  and  the  weighting  function 
peaks  at  the  600  mb  level.  Later  in  this  report  it  will  be  shown  how 
water  vapor  and  underlying  surface  effects  contribute  adversely  to  re- 
trieving atmospheric  temperature  profiles  fiom  satellite  radiometric  data. 
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3.1.2  L^rge  Data  Basi 


In  Older  to  St-'aicn  ■ ■:  i.  i.  , n.i.nii  :h  ittuLutnble  to 

hardware  features  of  thr-  iw"  d.;  ; - , . -v  ■ . i=  aji  le  to  use  a 

largo  data  base  to  ^.'nhwuuj-  ue.  - i..  ..  a : /eruj  sin.rlar 

types  of  atmosphere.s . D iV.  ,i  ii  .e  e.dings  ovr  i the 

northern  hernispnere  w..r<.'  u,-..  o : ,•1,.  1..  aulard  deviation 

of  radiances  measui\-d  in  ^a'-n  m •...n.m  Ls  .1  tiie  two  birds. 

These  data  were  separated  into  -!  ; ■ it  ina  t,  Kingitude  regions. 

Tne  actual  number  of  cases  in  • aoii  r ud  • r ■ i n m is  given  in  Table  3.2. 

In  general,  the  distribution  ot  I'ai'UL-  ...  ..:.,rii.  im-'-ng  trie  longitude  sub- 

divisions for  eacn  bird.  Tins  is  ..oi  L;.<  e.i  ... , :.owcvoi  , for  the  latitude 
subdivisions.  Theie  are  aiioul  t.v.ci  :i.an.  troj.icai  as  mici-latituae 
cases  and  about  twice  as  many  mid-latiu.'Un  .is  enl.ir  cases.  Over  the 
entu-e  northern  hemisn.ieia  , tiii.  t. .'t  il  ;.i.:’;...or  ; o.isos  arc  divided  nearly 
evenly  between  Birds  8531  aiu:  ■-.30.  Ai.  inaivsis  .it  ttiese  data  is  given 


later  in  section  3.2. 
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Table  3.2  Number  of  Cases  Used  to  Generate 
Mean  Radiances  Measured  as  a 
Function  of  Bird,  Latitude,  and 
Longitude,  Data  Period  14-28  February 
1975. 
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3.1.3  Coiiicid('nl  Sdtollitc'  .ind  RAOI^  D.Ua  Rase 


In  order  to  do  an  indepth  search  for  the  major  error  sources 
in  the  satellite  tempeiature  retrieval  process,  a search  was  made  through 
the  more  than  3000  satellite  sounding  cases  to  select  only  those  where 
coincident  RAOB's  would  be  available  for  verification  purposes.  Out  of 
this  large  data  sample,  only  10  cases  could  be  found  that  matched  our 
coincident  requirement,  i.e.  with  the  satellite  and  RAOB  sounding  within 
60  nautical  miles  and  one  hour  in  time  of  one  another. 

Table  3.3  gives  the  geographical  location,  bird  number,  and 
type  of  low  level  temperature  profile  for  the  10  coincident  soundings 
within  the  northern  hemisphere.  Each  bird  has  5 cases  of  nearly  simul- 
taneous satellite  and  RAOB  soundings.  There  is  one  case  in  the  tropics , 

4 cases  in  the  mid-latitudes,  and  5 cases  in  the  polar  regions.  All 
cases  are  over  the  ocean.  No  observations  were  available  for  the  North 
Atlantic  Ocean  but  a variety  covering  the  Pacific  and  Arctic  Ocean  and 
the  Mediterranean  Sea  was  obtained. 

Plots  of  the  observed  RAOB  temperature  profiles  are  shown  in 
Figs.  3.2  and  3.11,  inclusive,  for  the  10  soundings.  For  reference 
purposes,  the  AFGWC  12-hour  forecast  temperatures  at  mandatory  pressure 
levels  arc  included  since  these  arc  the  values  that  arc  used  as  the 
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Table  3.3  Location  and  Number  of  Research  Profiles  for  Three  Categories 
of  Low  Lev'el  Atmospheric  Stability 


Region 

Dr\’ 

Adiabatic 

South  of  Japan 

34. 5N 

138. 5E 

9532 

Profiles 

East  of  Japan 

38.  IN 

141 . 8E 

9532 

Uipse 

Near  Wake  Island 

19. 3N 

166.  IE 

8531 

Profiles 

Mediterranean  Sea 

34.  IN 

33.  IE 

8531 

Sea  of  Japan 
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9532 

Near  Finland 
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8531 
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Barents  Sea 
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Based 

Baffin  Bay 

69. 5N 

81 . OW 

9532 

Inversion 

North  of  Alaska 

71 .9N 

155.5W 

8531 

Profiles 

North  of  Russia 

79. 8N 

75. 4E 

8531 
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initial  guess  temperature  profile.  These  LO  cases  were  classified  according 

I 

to  the  low  level  atmospheric  stabilitv'  conditions  between  the  350  and 
1000  mb  pressure  levels.  Notice  that  there  are  four  cases  where  a very 
stable  low  level  atmosphere  exists  due  to  ground  based  temperature  in- 
version conditions.  Two  cases  exhibit  relatively  unstable  conditions 
associated  with  dry  adiabatic  typo  profiles.  The  remaining  four  cases 
conform  with  the  more  standard  type  atmospheres  found  with  lapse 
profiles . 
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figure  i.Z  Observed  und  Forecast 
Temperature  Profile  - 
68.2  N Latitude  53.  (i  f 
Longitude 

February  21  , 1975,  0014Z 


ii  I 


Pressure  (mb) 


Degrees  Celsius 

Figure  3.4  Observed  and  Forecast 
Temperature  Profile  - 
71 .9  N Latitude  155.5  W 
Longitude 

February  19,  1975,  2337Z 
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Figure  3.8  Observed  and  Forecast 
Temperature  Profile  - 
19 . 3 N Latitude  166.1  E 
Longitude 

February  16,  1975,  0118Z 
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Figure  3.10  Observed  and  Forecast 
Temperature  Profile  - 
47.4  N Latitude  141.7  E 
Longitude 

February  20,  1975,  2142Z 
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3.2  MEASURED  RADIANCES  EOR  BIRDS  8531  and  9532 


If  a large  enough  data  sample  exists,  it  should  be  possible  to 
analyze  the  measured  radiances  at  each  channel  from  one  satellite  and 
compare  it  with  those  from  another.  It  the  center  frequency  of  a partic- 
ular channel  was  off  or  had  drifted,  or  if  the  filters  had  changed  so 
that  the  band  width  allowed  a different  energy  spectrum  to  impinge  upon 
the  sensor,  or  if  other  satellite  hardware  rspects  changed,  then  it  would 
be  difficult  to  expect  that  such  changes  would  be  the  same  for  two 
different  satellites.  An  analysis  was  made  of  the  mean  and  standard 
deviation  of  the  radiances  detected  for  each  bird  as  a function  of  sensor 
channel,  latitude,  and  longitude.  In  particular,  comparisons  were  made 
to  determine  the  hardware  responses  for  each  sensor  channel  from  each 
bird.  Even  though  one  bird  is  a morning  and  another  is  a noon-time 
sensor,  a serious  deficiency  or  equipment  degradation  should  produce 
noticeable  differences  in  the  measured  radiances  for  a particular  channel 
when  equipment  discrepancies  exist  between  the  two  birds.  Furtherir.ore , 

before  the  adverse  effects  of  such  things  as  water  vapor  or  surface 
effects  can  be  placed  within  a proper  context,  it  is  important  to  know 
the  mean  conditions  expected  for  a particular  channel  and  the  normal 


variance  about  that  mean. 


A representative  measure  of  the  mean  radiance  and  standard 
deviation  as  a function  of  channel  and  latitude  is  given  in  Table  3.4 
for  Bird  8531.  In  general,  channel  3 which  has  a weighting  function 
which  peaks  closest  to  the  tropopauso  has  the  lowest  power  (or  tempera- 
ture) and  shares  with  channel  4 the  lowest  standard  deviation  of 
measured  radiance.  As  would  be  expected,  the  measured  radiant  power 
beginning  with  channel  4 or  greater  (with  weighting  functions  peaking  at 
400  mb  or  greater)  decreases  from  the  tropics  to  the  mid-latitudes  to 
the  polar  regions.  The  standard  deviation  of  satellite  measured  radiance 
within  the  tropics  is  about  half  of  that  for  extra-tropical  latitudes.  This 
should  be  expected  since  the  most  significant  day  to  day  atmospheric 
changes  take  place  within  the  extra-tropical  latitudes. 

Since  the  tropics  are  relatively  quiescent  and  do  not  play  an 
impoitant  roh>  in  most  day  to  day  numerical  weather  forecasts,  it 
oecomes  important  to  compare  changes  that  occur  between  mid-latitudes 
and  polar  regions.  Tor  the  more  than  1200  soundings  of  mid- latitudes 
and  polar  regions,  th(  variance  of  the  radiance,  measured  at  each 
channel  that  peaks  within  the  troposphere,  is  nearly  the  same.  That  is 
to  say,  that  whatever  diffi’renccs  exist  between  mid-latitude  and  polar 


TROPICS  - 0 to  2b  North  Latitude 


Channel 

J-  ^ A A ^ ^ L ^ 

Mean  Radiance  57.54  44.82  40.88  54. b5  76.31  92.33  117.8  123.6 

Standard  Deviation  1.47  1.69  .65  2.7  2 3.83  3.75  2.63  4.27 

MID-LATITUDE  - 26  to  62  North  Latitude 

Channel 

illlAlZ  8 

Mean  Radiance  56.38  47.83  44.27  49.8  63.81  75.67  91.58  105.22 

Standard  Deviation  5.  5.06  3.63  2.92  5.8  7.92  13.18  8.06 

POLAR  - 62  to  9 0 Nortn  Latitude 

Channel 

iZiilZIB 

Mean  Radiance  45.5  37.95  35.55  40.65  50.76  58.52  61.56  87.27 

Standard  Deviation  4.97  5.42  3.87  2.74  4.29  6,93  13.27  8.99 

NOR'n-IERN  HEMISPHERE  - 0 to  90  North  Latituae 

Channel 

12^456  I 8 

Mean  Radiance  55.59  44.65  41.02  51.47  69.58  83.39  103.32  114.14 

Standard  Deviation  5.25  4.6  3,54  5.53  10.24  13.42  21.8  14.44 


Table  3.  ! Measured  Radiance  Mean  and  Standard 
Deviation  Values  for  Bird  8531  for  0 to 
180  Degrees  East  Longitude  as  a Function 
of  Latitude  and  Sensor  Channel 


7') 


regions,  satellite  radiometers  detect  the  same  signal  power  variations. 

This  fact  will  become  increasingly  more  imporant  when  the  effects  due 
to  watei  vapor  contributions  are  discussed  in  section  3.4,  where  it  will 
be  shown  that  water  vapor  cannot  possibly  account  for  the  variance  of 
radiance  observed  in  the  polar  regions. 

When  considering  all  latitudes  and  longitudes,  it  can  be  shown, 
as  in  Fig.  3.12,  that  the  average  measured  radiance  for  each  Bird  is 
highly  correlated  with  the  other,  with  a correlation  coefficient  of  0.9978. 
This  im.plies  that  within  the  overall  gross  sense,  observations  from  one 
satellite  correspond  closely  with  the  other  satellite. 

A comparison  was  made  betwet  n the  average  measured  difference 
in  radiance  between  Birds  8531  minus  9532  as  a function  of  sensor 
channel  and  latitude  belt,  as  shown  in  Figs.  3.13,  3.14,  and  3,15.  In 
all  cases,  for  all  latitudes,  channel  4 radiance  for  Bird  8531  was  lower 
(colder)  than  that  for  Bird  9 532.  This  feature  proved  to  be  the  exception, 
ho'.vever.  As  latitude  of  observation  increased,  so  did  the  difference 
between  measured  radiance  between  Birds  8531  and  9532.  In  general. 

Bird  8531  received  more  radiant  power,  showed  warmer  temperatures,  than 
Bird  9532  for  ever  increasing  latitudes.  In  fact.  Fig.  3.16  shows  a 
composite  of  all  satellite  soundings  at  all  latitudes  and  longitudes  for 
each  bird  as  a function  of  sensor  channel.  Once  again,  channel  4 is 


no 


odd,  showing  that  the  sensor  for  Bird  8531  reports  lower  or  colder 
radiances  than  the  same  sensor  for  Bird  9532.  On  tne  other  hand,  all 
remaining  sensor  channels  of  Bird  85  31  report  higher  or  warmer  values 
than  those  from  Bird  9532.  These  results  suggest  that  equipment  differ- 
ences or  degradations  need  to  be  explored  between  satellite  Birds  8531 
and  9532. 
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Figure  3.12  Correlation  of  Radiance 

Measured  from  Birds  8531 
and  9532 
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Figure  3.13  Difference  in  Radiance 

Measured  by  Both  Satellites 
in  Tropics 
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Figure  3.  14  Difference  in  Radiance 

Meajured  by  Both  Satellites 
in  Mid-Latitudes 
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Figure  3.16  Difference  in  Radiance 

Measured  by  Both  Satellites 
in  Northern  Hemisphere 
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3.3  DISCRIMINATION  CRITERIA 


Several  types  of  control  parameters  are  used  in  the  AECWC 
computer  program  SOUNDER  that  affect  acceptance  or  rejection  of  a 
satellite  sounding.  These  control  parameters  were  studied  to  determine 
their  influence  on  the  number,  type,  and  quality  of  the  satellite 
radiance  observations  that  are  accepted  for  use  in  the  inversion  process 
Section  2.1.3,  the  subprogram  ITERAT,  discusses  in  detail  functions  of 
the  control  parameters  DELVVND,  CLDDIF,  and  HOTSFC.  DELWND  is 
the  radiance  measured  for  window  channel  7 minus  the  radiance  cal- 
culated from  both  the  underlying  surface  and  atmosphere.  The  1000  mb 
air  temperature  is  used  as  the  "ground"  temperature  to  compute  surface 
radiation  contriliutions . The  value  of  CLDDIF  is  used  to  test  for  clear 
radiance  and  to  reject  a sounding  where  clouds  might  be  causing  a 
difference  between  observed  and  computed  radiance  values.  The  value 
of  HOTSFC  is  used  to  reject  a sounding  when  a hot  surface  exists  or 
a bad  first  guess  temperature  profile  is  made. 

Analyses  were  made  on  the  same  data  base  using  values  of 
these  parameters  that  wore  operational  in  February  1975  and  that  are 
operational  today.  Results  show  that  SOUNDER  with  the  present  values 
accepts  only  1/4  as  many  cases  for  completing  an  inferred  profile  as 


it  did  in  197  5.  This  was  caused  by  imposing  a more  restrictive  accep- 
tance criteria  of  K1  tor  both  CLDDIT  and  HOTSFC  in  the  hopes  of 
reducing  rattier  large  retrieval  errors. 

Let  us  consider  the  window  channel  difference  values  of 
DELWND  in  Table  3.5  for  the  10  research  profiles  having  coincident 
satellite  and  RAOB  soundings.  In  order  for  any  of  these  soundings  to 
be  acceptable  for  use  in  the  temperature  retrieval  process,  both  the 
low  level  cloud,  CLDDIF,  and  hot  surface,  HOTSFC,  criteria  would 
have  to  equal  or  exceed  these  DELWND  values.  With  today's  opera- 
tional criteria  of  1.1,  there  is  only  one  case  in  ten  that  would  be 
accepted  cut  of  the  sample  where  the  nearest  (coincident)  RAOB  was 
used  to  calculate  surface  and  atmospheric  radiance  in  the  window 
channel.  Furthermore,  that  one  case  exists  in  the  tropics  where 
temperature  profile  data  are  usually  not  as  critical  to  the  success  of 
numerical  analysis  and  prediction  models.  Using  the  12-hour  forecast 
data  in  computing  the  radiances  to  obtain  DELWND,  Table  3.5  values 
show  only  three  of  the  ten  cases  would  be  acceptable  to  the  re- 
trieval process.  All  three  of  the  acceptable  cases,  however,  have 
relatively  standard  type  lapse  profiles.  It  appears,  therefore,  that  the 
discrimination  criteria  in  use  today  are  eliminating  those  cases  having 
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Tabic  3.5  Low  Level  Clouds  and  Hot  Surface  Discrimination 
Values  for  the  Ten  Research  Profiles 
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the  most  significant  types  of  meteorological  profiles. 


Striking  results  are  shown  in  Fig.  3.17  where  DELWTMD  values 
are  plotted  against  the  observed  RAOB  temperature  gradient  from  1000 
to  850  mb  level  for  each  of  the  ten  research  cases.  Notice  first  that 
the  distribution  of  plots  is  quite  different  for  Bird  9532  versus  8531. 
This  IS  duo  in  part  to  that  fact  that  Bird  8531  has  more  lapse  profile 
type  cases.  But  of  most  importance  is  the  17  mW  m”^  Sr~^  cm^  range 
in  radiance  associated  with  the  type  and  intensity  of  low  level  atmos- 
pheric stability.  This  range  m radiance  is  highly  significant  especially 
since,  as  shown  in  the  next  section,  it  greatly  exceeds  the  range  of 
radiances  associated  with  water  vapor  contributions  in  mid  and  high 
latitudes  and  is  comparable  with  the  maximum  extremes  possible  within 
tlie  tropics.  It  IS  not  the  low  level  temperature  gradients  per  se  but 
rather  the  effects  due  to  the  underlying  surface  that  are  important. 

These  surface  effects  are,  however,  portrayed  by  the  low  level  atmos- 
pheric stability  conditions. 

When  a ground  based  temperature  inversion  exists,  the  air 
temperature  increases  with  height.  Also  under  these  conditions,  the 
underlying  surface  is  colder  than  the  air  temperature.  The  reverse  is 
true  when  the  atmosphere  becomes  more  unstable  and  the  low  level 
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Figure  3.17  Discrimincition  Criterion  DELWND 
Related  to  Low  Level  Atmospheric 
Stability 
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temperature  profile  approaches  or  exceeds  dry  adiabatic  conditions. 

Since  the  surface  radiance  contributions  are  computed  using  the  1000  mb 
level  air  temperatures,  significant  differences  can  be  expected  between 
measured  and  calculated  when  the  1000  mb  level  air  temperature  is  non- 
representative  of  the  surtoice  temperature.  This  becomes  more  pronounced 
as  the  low  level  atmospheric  stability  deviates  from  the  standard  lapse 
type  profiles.  Thus  it  can  be  seen  in  Tig.  3.17  that  for  ground  based 
inversions,  stable'  atmospheres  with  a negative  change  in  temperature 
with  pressure,  the  measured  satellite  radiance  in  the  window  channel 
IS  lower  (colder)  than  the  calculated  ground  and  atmospheric  contribu- 
tions giving  a negative  value  of  DELWKD.  The  reverse  is  true  for 
those  cases  wnen  the  atmosphere  becomes  increasingly  more  unstable. 
These  striking  effects  due  to  the  underlying  surface  are  discussed  in 
greater  detail  in  section  3.5. 

Tor  our  purpose  here,  it  must  be  stated  that  the  existing 
operational  values  for  CLDDIF  and  HOTSFC  are  too  restrictive.  For 
example,  when  the  measured  surface  radiance  is  lower  (colder)  than 
that  calculated  by  an  amount  exceeding  CLDDIF,  then  no  temperature 
retrievals  are  presently  made  from  satellite  data.  The  rationale  being, 
if  a very  low  level  cloud  exists,  which  could  not  be  detected  via 


previous  discrimination  procedures,  then  the  cloud  top  would  produce 
a lower  measured  radiance  in  the  window  channel  than  is  calculated 
and  thus  a cloud -free  line  of  sight  does  not  exist  and  a temperature 
retrieval  should  not  be  attempted.  As  shown  in  Fig.  3.17,  this 
rationale  eliminates  all  the  stable  cases  with  low  level  ground  based 
inversions  where  the  underlying  surface  air  temperature  can  be  con- 
siderably lower  (colder)  than  the  1000  mb  air  temperature.  In  fact, 
under  these  conditions,  if  a low  level  cloud  did  exist,  the  cloud  top 
temperature  would  be  warmer,  not  colder,  tnan  the  underlying  surface. 
For  several  research  profiles  with  ground  .lased  inversions,  the  ground 
air  temperature  is  colder  than  any  ambient  air  temperature  up  to  10,000 
feet  (700  mb).  Thus  restrictive  values  for  CLDDIF  eliminate  tne  possi- 
bility of  temperature  retrievals  under  these  conditions.  There  are 
times  when  CLDDIF  criteria  could  cTi.minate  low  level  cloud  cases  but 
tliese  times  can  only  exist  during  conditions  of  lapse  type  profiles. 
Since  the  first  guess  temperature  profile,  given  by  the  12-hour  fore- 
cast field,  provides  a measured  of  the  low  level  atmospheric  stability, 
It  could  be  used  to  improve  the  discrimination  or  selection  process  for 
satellite  soundings  worthy  of  temperature  retrievals. 

Table  3.6  shows  the  different  combinations  of  criteria  for 


CLDDIF,  HOTSFC,  SNOIS,  and  EPS  used  within  the  GAC  research 


SOUNDER,  Frr  in  Table  3.5  it  can  be  seen  that  a value  of  5 for  CLDDIF 

and  HOTSFC  will  allow  all  cases  for  Bird  8531  to  be  retrieved.  A larger 

value  was  necessary  to  accomplish  temperature  retrievals  for  all  Bird 

9532  cases.  The  purpose  here  was  to  insure  that  all  10  research  profile 

cases  were  acceptable  for  temperature  retrievals  in  order  to  determine 

the  effects  of  SNOIS  and  EPS  on  the  resultant  inverted  temperatures. 

SNOIS  is  the  ratio  of  signal  to  noise  of  the  instrument  and  is  used  to 

determine  if  the  sum  of  the  square  of  the  deviations  between  measured 

and  calculated  radiances  is  small  enough  to  allow  convergence  for  a 

solution.  EPS  is  the  epsilon  value  used  to  monitor  the  maximum  difference 

allowable  between  calculated  and  measured  radiance  from  any  channel. 

Table  3.6  Cloud  and  Hot  Surface  Discriminants  and  Solution 
Convergence  Values 

For  Bird  4953  2 (5  locations) 


CLDDIF  5 

HOTSFC  5 

SNOIS  3 

EPS  1 


5 20 

5 20 

1.5  3 

0.5  1 


20  20 

20  20 

0.3  0.1 

0.1  0.033 


For  Bird  #853  1 (5  locations) 


CLDDIF  5 5 5 

HOTSFC  5 5 5 

SNOIS  1.5  0.3  0.1 

EPS  0.5  O.l  0.033 


These  control  parameters  and  how  they  interface  with  SOUNDER  are 
discussed  in  detail  in  section  2.1.3.  In  general,  it  can  be  said  that 
minor  differences  were  obtained  in  the  retrieved  temperature  profiles 
even  though  a large  range  of  values  was  used  for  SNOIS  and  EPS.  Even 
for  a five-fold  difference  in  SNOIS  and  EPS  values,  as  shown  in 
Fig.  3.18,  temperature  retrievals  at  different  pressure  levels  differed  by 
only  a few  tenths  of  a degree.  On  the  average,  the  largest  differences 
occurred  in  the  tropopause  region  and  at  about  the  70  mb  level.  None 
of  these  differences  could  account  for  the  source  of  the  larger  errors 
in  the  satellite  retrieval  process. 
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Figure  3.18  Average  Difference  in  Retrieved 

Temperature  for  Two  Sets  of  SNOIS  and 
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SNOIS  = 1.5  and  EPS  = 0.5) 
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3.4  WATER  VAPOR  CONTRIBUTIONS 


No  attempt  is  made  in  the  existing  operational  system  to  in- 
clude variable  effects  due  to  water  vapor.  Recall  from  Table  2.5  that 
a fixed  relative  humidity  value  is  assigned  to  the  28  levels  extending 
from  the  surface  to  the  200  mb  level.  Typical  assigned  relative  humidity 
values  are  about  73%  near  the  surface,  65%  at  850  mb,  55%  at  700  mb, 
47%  at  500  mb,  32%  at  300  mb,  and  16%  at  200  mb  level.  For  fixed 
relative  humidity  values,  as  the  first  guess  temperatures  increase  or 
decrease,  so  does  the  water  vapor  mixing  ratio,  which  is  so  important 
in  determining  the  concentration  of  water  vapor  molecules  which  con- 
tribute to  the  final  measured  or  calculated  radiance.  Via  this  feature, 
the  warmer  low  level  temperatures  in  the  tropics  account  for  much 
larger  water  vapor  contributions  than  do  the  colder  polar  temperatures. 

In  this  manner  the  existing  SOUNDER  routine  accounts  for  major  latitude 
variations  of  water  vapor  effects. 

A recent  study  by  Weinreb  investigated  the  effects  of  water 
vapor  contributions  on  the  retrieved  temperature  profile.  It  was  shown 
that  "changes  in  the  assumed  mixing  ratio  affect  the  retrieved  tempera- 
ture chiefly  through  changes  in  the  calculated  radiances".  Furthermore, 
it  was  shown  that  when  the  initial  forecast  guess  of  the  temperature 
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profile  was  close  to  the  actual  that  errors  in  water  vapor  estimates 

resulted  in  larger  errors  in  the  minimum-information  technique  than  the 

initial  forecast  guess.  On  the  other  hand,  when  the  initial  forecast 

temperature  profile  guess  had  large  errors,  the  minimum-information 

technique,  even  with  water  vapor  errors,  produced  temperature  profiles 

that  were  an  improvement  upon  the  initial  forecast  guess.  In  essence, 

Weinreb  suggests  that,  for  the  most  part,  errors  in  estimating  water 

vapor  contributions  are  principal  factors  in  attaining  accurate  satellite 

temperature  retrievals.  In  order  to  place  these  statements  in  perspective 

with  the  underlying  surface  effects  discussed  in  the  next  section,  the 

changes  in  the  computed  radiances  need  to  be  considered  for  errors  in 

the  water  vapor  mixing  ratio.  Weinreb' s computational  results  from  a 

Midway  Island  RAOB  shows  a linear  relation  between  radiance  changes 

and  errors  in  the  mixing  ratio.  A 50%  error  in  mixing  ratio  produced 

-1  -2  1 

changes  in  radiance  of  3.6,  2.6,  1.1,  and  0.2  mW  sr  m cm  for 
sensor  channel  7 (835  cm  ^),  6 (747  cm  ^),  5 (725  cm  ),  and  4 
(708  cm"^),  respectively.  It  will  be  shown  that  these  radiance  changes 
are  small  in  comparison  with  underlying  surface  effects. 

The  maximum  effects  possible  due  to  water  vapor  on  the 
satellite  temperature  retrievals  were  studied.  Radiances  were  calculated 
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for  each  of  the  10  research  l^OB  profiles  assuming  complete  saturation 


(wet)  at  all  levels  in  the  atmosphere.  Tltcse  computations  were  repeated 
assuming  zero  moisture  (dry)  at  all  levels.  Figure  3.19  presents  the 
difference  in  the  temperatures  retrieved  via  the  minimum  information 
technique  using  totally  wet  and  dry  conditions  imposed  upon  the  same 
RAOB  sounding.  The  largest  temperature  retrieval  errors  due  to  water 
vapor  occur  in  the  850  and  7 00  mb  regions.  An  oscillation  exists  in 
the  vertical  distribution  of  retrieved  temperature  differences  such  that 
the  surface  value  is  nearly  the  same  sign  and  magnitude  as  at  100  mb. 
Whereas  the  maximum  positive  differences  around  7 00  mb  give  way  to 
negative  differences  about  the  300  mb  level.  Thus  it  can  be  seen  how 
errors  in  the  lower  portion  of  the  atmosphere  propagate  to  higher  levels. 

Very  large  differences  (22.5°C)  exist  in  temperatures  retrieved 
at  700  mb  from  a saturated  and  totally  dry  tropical  sounding,  as  shown 
in  Fig.  3.19.  Mid-latitude  differences  are  about  half  of  the  extrem.e 
values  found  in  the  tropics.  Even  in  the  polar  regions  where  cold 
temperatures  inhibit  high  water  vapor  mixing  ratios,  a 2°C  difference 
can  be  expected  between  retrieved  temperatures  from  a wet  and  dry 
profile. 

The  extreme  effects  of  a totally  saturated  versus  a completely 


Pressure  (mb) 


Retrieved  Temperature  Difference  from 
Wet  minus  Dry  Profile 

Figure  3.19  Levels  in  the  Atmosphere  Where  the 

Retrieved  Temperatures  are  Influenced  by 
Water  Vapor 


1(X) 


dry  atmosphere  on  calculatea  radiance  are  shown  in  Fig.  3.20.  Radi- 
ances calculated  from  a dry  minus  those  from  a wet  profile  are  plotted 
against  the  average  of  the  saturation  mixing  ratios  at  mandatory  levels 
from  1000  mb  to  500  mb  level.  The  radiance  difference  for  channels  4, 
5,  6,  and  7 are  plotted  at  three  average  saturation  mixing  ratios  that 
were  derived  from  a polar,  mid-latitude  and  tropical  RAOB  profile. 

These  radiance  differences  arc  the  extremes  possible  for  water  vapor 
effects.  The  greatest  radiance  difference  (19  mW  m sr  cm)  is 
found  in  window  channel  7 for  a tropical  profile.  In  general,  the 
radiance  difference  due  to  moisture  extremes  decreases  as  the  channel 
numoer  decreases  and  the  depth  of  atmosphere  under  interrogation  de- 
creases.  Radiance  differences  of  less  than  5 mW  m ^ sr  cm  occur 
at  mid-latitudes  for  all  channels.  These  values  can  produce  pronounced 
errors  in  retrieved  temperatures  and  efforts  should  be  made  to  incorpo- 
rate satellite  derived  moisture  information  into  the  retrieval  process. 
Moisture  effects  arc  nearly  non-cxistent  for  the  polar  regions  where 

the  average  saturation  mixing  ratio  between  1000  and  500  mb  is  less 

-7  - 1 

than  one  and  the  radiance  difference  is  less  than  0.4  mW  m sr  cm. 
These  radiance  differences  from  a dry  and  wet  profile  will  be  shown 
to  be  small  for  mid-  and  high-latitudes  when  compared  with  underlying 
surface  effects. 
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Between  1000  and  500  mb  Level 
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3.5  UNDERLYING  SURFACE  ErEECTS 


Section  3.3  described  the  physics  and  rationale  for  effects  of 
underlying  surfaces  on  calculated  radiances.  In  short,  it  says  that  for 
stable  atmospheres  with  ground-based  inversions,  air  temperature  in- 
creases with  height  so  that  the  surface  temperatures  are  colder  than 
the  1000  mb  level  air  temperature.  The  reverse  is  true  for  the  unstable 
dry  adiabatic  type  conditions.  Presently  the  1000  mb  air  temperature 
is  used  to  compute  surface  contributions  to  satellite  radiances.  But 
due  to  the  nonrepresentativeness  of  these  air  temperatures  in  depicting 
surface  temperatures,  the  largest  single  source  of  error  in  calculated 
radiances  is  due  to  these  underlying  surface  effects. 

The  ten  research  RAOB  profiles  were  separated  into  three 
categories  of  ground-based  inversions,  dry  adiabatic,  and  the  more 
standard  lapse  type  profiles.  In  order  to  determine  surface  effects  on 
each  sensor  channel,  the  expected  radiance  was  calculated  using  first, 
the  coincident  RAOB  and  then  the  12-hour  forecast  temperature  data. 
Differences  were  obtained  between  measured  and  calculated  radiances 
and  grouped  according  to  the  three  stability  categories.  These  radiance 
differences  were  added  separately  for  each  category  to  obtain  the 
algebraic  and  the  absolute  mean  differences,  which  are  tabulated  in 


r 


USING 

NEAREST 

RAOB 

DATA 


USING 

TWELVE  HOUR 
FOREGAST 
DATA 


Average 

Absolute 

Average 

Absolute 

Mean 

Mean 

Mean 

Mean 

Difference 

Difference 

Difference 

Difference 

DRY 

1 .214 

.216 

.384 

.384 

ADIABATIG 

2 1.706 

1.706 

4.209 

4.209 

PROFILES 

0) 

c 

3 .541 

.541 

2.102 

2.102 

4 .163 

.289 

2.231 

2.231 

fO 

5 3.098 

3.098 

2.732 

2.732 

O 

6 3.546 

3.546 

.119 

.832  i 

7 8.710 

8.710 

-1.488 

1.488 

LAPSE 

1 - .209 

1.844 

.582 

1.699 

PROFILES 

2 .665 

2.079 

1.585 

2.344  ! 

'o 

3 .162 

1.016 

.804 

.804 

c 

c 

4 - .740 

1 .154 

- .507 

1.004  ; 

fO 

x: 

5 1.645 

1.645 

2. 109 

2.109 

O 

6 .310 

.804 

.989 

.989  ; 

7 .626 

2.752 

1.766 

2.304  i 

GROUND 

1 

-2.865 

2.865 

-2.575 

2.575 

BASED 

2 

-1.393 

1.393 

- .882 

1.593 

INVERSION 

a 

c 

c 

3 

-1.558 

1.568 

- .372 

1.643 

PROFILES 

4 

-1.864 

1.864 

- .443 

.974 

r* 

5 

-1.528 

1.689 

.056 

.623 

o 

6 

-2.540 

2.540 

- .512 

1.249 

7 

-3.799 

3.799 

-1.045 

2.470 

Table  3.7  Difference  between  measured  and  calculated  radiances  for 
three  types  of  low  level  temperature  profiles.  Galculated 
radiances  were  made  separately  using  nearest  RAOB  data 
and  twelve-hour  forecast  data. 
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Table  3.7,  Thus  if  the  average  mean  difference  has  exactly  the  same 
magnitude  as  the  absolute  moan  difference,  then  differences  for  each 
individual  case  within  the  ensemble  must  have  the  same  sign.  That  is 
to  say  that  all  cases  are  influenced  the  same  way.  This  feature  is 
exemplified  for  almost  every  sensor  channel  for  both  ground-based  in- 
version and  dry  adiabatic  profiles  when  RAOB  data  are  used  to  compute 
radiances.  Notice  for  ground-based  inversion  conditions,  negative 
values  of  the  average  mean  difference  show  that  satellite  measured 
radiance  is  lower  (colder)  than  that  calculated.  Furthermore,  and  this 
is  very  significant,  these  effects  arc  prominent  in  all  seven  channels, 
including  channel  1 (668.5  cm  on  the  peak  of  the  CO2  band  where 
the  weighting  function  supposedly  terminates  well  above  the  under- 
lying surface.  The  channel  1 average  mean  difference  of  -2.865 
mW  m sr  cm  for  ground-based  inversion  cases  in  the  polar  regions 
is  more  than  half  the  observed  standard  deviation  value  that  represents 
the  variability  for  all  518  polar  cases. 

Influence  of  the  underlying  surface  is  even  nore  pronounced 
for  channels  5,  6,  and  7 for  dry  adiabatic  cases  which  occurred  in 
mid-latitudes.  All  seven  channels  show  the  measured  radiance  to  be 
higher  (warmer)  than  that  calculated.  This  is  due  to  the  underlying 


surface  being  hotter  than  the  1000  mb  air  temperature.  An  average  mean 


difference  between  measured  and  RAOB  calculated  radiances  for  channel  7 

-2-1  9 

is  8.710  mW  m sr  cm  for  dry  adiabatic  cases  and  -3.799  mW  m“‘^ 

sr~^  cm  for  ground  based  inversion  cases.  This  gives  a range  in  the 
mean  values  of  12.509  mW  m”^  sr“^  cm.  To  place  this  range  value  in 
perspective,  it  can  be  seen  from  Fig,  3.20,  that  this  value  due  to 
underlying  surface  effects  is  more  than  twice  as  large  as  that  attainable 
from  extreme  mid-latitude  water  vapor  effects,  i.e.  totally  saturated 
versus  a totally  dry  profile.  An  analysis  for  channels  5 and  6 shows 
similar  results.  Thus,  it  can  be  concluded  that  the  underlying  surface 
effects  are  more  important  than  water  vapor  effects  in  interpreting  in- 
formation from  radiance  measurements  from  channels  5,  6,  and  7. 

Properly  accounting  for  these  underlying  surface  effects  should  greatly 
improve  the  ability  to  remotely  measure  low  level  temperature  profiles. 

Even  for  standard  lapse  type  profiles.  Table  3.7  shows  measured 
radiances  on  the  average  are  higher  (warmer)  than  those  calculated.  Al- 
though surface  effects  are  recognizable  and  consistent  for  lapse  profiles 
they  are  not  as  striking  as  those  for  ground-based  inversion  or  dry 
adiabatic  conditions.  Results  using  12-hour  forecast  temperature  data 
to  compute  radiances  show  characteristics  similar  but  not  as  pronounced 


1(K. 


as  those  using  RAOB  data. 
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Low  Level  Temperature  Gradient  (AT/AP)(°C/mb)(RAOB) 

Figure  3.21  Radiance  Differences  as  a 
Function  of  Atmospheric 
Stability  - Channel  One 


Channel  2 


Low  Level  Temperature  Gradient  (AT/AP)(°C/mb)(RAOB) 

Figure  3.22  Radiance  Differences  as  a 
Function  of  Atmospheric 
Stability  - Channel  Two 
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Channel  4 


Low  Level  Temperature  Gradient  (AT/AP)(°C/mb)(RAOB) 

Figure  3.24  Radiance  Differences  as  a 
Function  of  Atmospheric 
Stability  - Channel  Four 
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Low  Level  Temperature  Gradient  (ATMP)(‘’C/mb)(RAOB) 

Figure  3.25  Radiance  Differences  as  a 
Function  of  Atmospheric 
Stability  - Channel  Five 


Measured  - Calculated  (RAOB)  Radiance 
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Low  Level  Temperature  Gradient  (AT/AP)('’C/mb)(RAOB) 

Figure  3.26  Radiance  Differences  as  a 
Function  of  Atmospheric 
Stability  - Channel  Six 
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Channel  7 
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Low  Level  Temperature  Gradient  ( T P)(^C  mb)(RAOB) 

Figure  3.Z7  Radiance  Differences  as  a 
Function  of  Atmospheric 
Stability  - Channel  Seven 
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Plots  were  made  to  show  the  relationship  between  measured 
minus  calculated  radiances  to  corresponding  individual  values  of  low 
level  atmospheric  stability  for  each  of  the  seven  sensor  channels,  see 
Figs.  3.21  to  3.27.  All  calculated  radiances  and  temperature  gradients 
from  1000  to  850  mb  were  made  from  RAOB  data.  Successively  viewing 
the  figures  illustrates  an  increasing  influence  of  surface  effects  on 
information  contained  with  increasing  sensor  channel  number.  A relation- 
ship even  exists  for  both  Bird  s channel  1 sensor  whose  weighting 
function  peaks  at  30  mb  level.  Beginning  with  channels  5,  6,  and  7 
where  the  weighting  functions  peak  at  7 00  mb,  1000  mb,  and  surface 
levels,  respectively,  the  difference  between  measured  and  calculated 
radiance  becomes  highly  correlated  with  low  level  temperature  gradients. 

The  window  channel  7 results  in  Fig.  3.27  show  a range  of  more  than 
-2  “1 

17  mW  m sr  cm  in  the  difference  between  measured  and  calculated 
radiances.  This  error  due  to  underlying  surface  effects  greatly  exceeds 
any  other  error  source  investigated  for  cloud  free  conditions. 

Underlying  surface  temperatures  contribute  to  outgoing  radia- 
tion more  strongly  than  accounted  for  in  operational  practice  today. 

More  importantly,  the  sign  and  magnitude  of  the  low  level  temperature 
profile  appear  suitable  as  aids  to  correct  for  these  deficiencies.  This 
IS  shown  graphically  in  Fig.  3.28  where  satellite  retrieved  temperature 
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Retrieved  Temperature  Errors  at  1000  mb  Level 


Low  Level  Temperature  Gradient  (AT/AP)(°C/mb) 

Figure  3.28  Retrieved  Temperature  Errors  at  1000  mb 
Level  Related  to  Forecast  Low  Level  Air 
Temperature  Gradients 
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errors  at  1000  mb  level  are  plotted  relative  to  corresponding  tempera-  | 

ture  gradients  provided  by  the  AFGWC  12-hour  forecast  for  850  and  ' 

1000  mb  levels.  Retrieved  temperatures  experienced  a 15°C  range  in 
errors  for  the  ten  research  profile  cases.  It  can  be  seen  that  retrieved 
temperature  errors  at  1000  mb  are  highly  related  (correlation  coefficient 
of  0.912)  to  forecast  temperature  gradients.  In  conclusion,  objective 
techniques  can  and  should  be  developed  to  account  for  large  under- 
lying surface  effects  on  an  operational  basis  at  AFGWC. 
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4 STATISTICAL  ANALYSIS 


4.1  INTRODUCTION 

The  stepwise  regression  or  subset  regression  analysis  was 
accomplished  at  Scientific  Systems,  Inc,,  under  subcontract  with  the 
Geo-Atmospherics  Corporation, 

4.2  OBJECTIVES 

The  objectives  of  the  project  which  relate  to  the  statistical 
analysis  problem  were: 

1,  Compare  operational  soundings  of  vertical  temperature  profiles  ob- 
tained from  the  infrared  Special  Sensor  (SSE)  package  aboard  the  Defense 
Meteorological  Satellite  Program  (DMSP)  with  coincident  radiosonde 
observations  and  available  forecast  and  analysis  field  data, 

2,  Investigate  the  cause  and  effect  of  errors  in  the  operational  soundings, 

3,  Develop  techniques  which  mdnimize  the  vertical  temperature  profiling  errors, 

4,  Perform  similar  analysis  for  computation  of  thickness  between  mandatory 
levels , 

4.3  STATISTICAL  DATA  BASE 
4,3,1  Data  Base  Structure 

The  data  base  selected  for  use  in  this  statistical  study  was 
obtained  from  the  data  available  at  the  Air  Force  Geophysics  Laboratory 
(AFGL) , The  data  covered  the  Northern  Hemisphere  during  a two-week  period 
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from  February  lA  to  28,  1975  and  contained  DMSP  first  guess  profiles, 

DMSP  final  temperature  estimates.  Navy  TTPACK  data,  and  RAOB  (radiosonde)  data. 
All  data  had  passed  the  tests  used  in  the  DMSP  SOUNDER  program  and  thus 
were  judged  as  "clear  column"  radiances. 

The  data  on  tape  represented  only  4%  of  the  measured  data. 

This  saved  data  is  the  data  for  which  there  is  the  greatest  possibility 
of  no-cloud  conditions,  so  that  the  dear-column  radiance  assumption  used 
in  the  minimum-information  solution  is  made  as  valid  as  possible.  Each 
measurement  consists  of  5 individual  radiance  measurements  per  scan  line 
with  5 scans  being  made.  One  of  these  25  measurements,  which  is  made  over 
the  clearest  field  of  view,  according  to  the  channel  4/channel  7 radiance 
ratio  test,  is  kept  for  use  in  the  temperature  or  height  retrievals. 

Statistical  data  compiled  at  AFGL  on  this  data  base  has  been 
analyzed  to  determine  the  characteristics  of  the  observed  errors.  In  com- 
paring TTPACK  vs.  SSE  inversions,  the  following  points  should  be  made. 

Large  errors  occur  predominantly  at  three  pressure  levels:  (1)  surface, 

(2)  tropopause,  (3)  upper  stratosphere.  The  errors  within  the  stratosphere 
(primarily  at  and  below  30  mb)  are  primarily  negative  and  appear  to  be 
correlated  with  an  isothermal  first  guess  profile  at  these  levels.  Errors  at 
the  tropopause  do  not  appear  to  be  related  to  those  in  the  upper  atmosphere. 
Tropopause  errors  appear  to  be  most  frequent  at  mid-latitude  soundings.  The 
mean  error  at  all  levels  tracks  the  inversion  coefficient  for  the  667/cm  ^ 
channel . 

In  order  to  generate  a data  base,  several  different  selection 
criteria  were  used.  In  order  to  minimize  the  effect  of  spatial  horizontal 
and  temporal  inhomogeneities,  it  is  desirable  that  the  DMSP  soundings  corres- 
pond closely  in  space  and  time  with  the  TTPACK  data.  Using  a threshold  of 
60  nautical  miles  and  3 hours,  a total  of  1067  matches  were  found  (563  for 
the  first  week  and  504  for  the  second  week).  A subsequent  set  of  matches 
were  generated  with  tighter  thresholds  but  the  number  of  matches  was  signifi- 
cantly reduced.  It  was  decided  that  the  1067  matches  were  large  enough  in 
number  to  make  a statistically  significant  test.  At  the  same  time,  the 
thresholds  were  felt  to  be  sufficiently  tight  to  minimize  the  effects  of 
spatial  and  temporal  inhomogenelties.  Therefore,  this  basic  data  base  was 
adopted  for  use  in  the  study. 


One  objective  of  this  study  was  to  assess  tiie  effect  of  latitude 
and  zenith  angle  on  t lie  regression  models.  Variations  in  latitude  strongly 
affect  the  temperature  profile,  while  variations  in  zenith  angle  jiroduce 
changes  in  the  observed  spectral  ratliances,  even  for  invariant  temperature 
prof iles. 

The  latitudes  were  partitioned  into  two  sets,  one  from  0 to 
45“N  and  one  from  45°N  to  90°N,  and  the  zenith  angles  (0)  were  also  parti- 
tioned into  two  sets,  one  for  0 < G<0.7  radian  and  one  for  0 > 0.7  radian. 

Only  two  latitude  regions  were  chosen  to  maintain  as  much  statistical  signifi- 
cance as  possible.  The  zenith  angle  partition  was  selected  to  give  sets 
of  approximately  equal  size.  It  is  curious  that  the  data  accepted  in  SOUNDER 
had  a relatively  large  number  of  high  zenith  angles.  This  is  suggestive  of 
an  unwanted  bias  in  tlie  computations  involved  in  making  the  cloud/no  cloud 


decision . 


Tlie  number  ol  matches  for  the  various  partitions  are  given  in 


Table  4.1. 


TABLE  4.1 


Number  of  Matches  between  DMSP  and  TTPACK 
for  Latitude  and  Zenith  Angle  Partitions 


He  sc  r i t>t  i on 


Week  1 

14-21  Feb.  1975 


Week  2 

22-28  Feb.  1975 


all  data  563  504 

0 < A < 45°N  . 164  332 

45°N  < A 9n°N  399  172 

0 < 0 < 0.7  rad.  303  275 

0 > 0.7  rad.  260  229 

4.3.2  Validation  of  Data  Base 

In  order  to  perform  a meaningful  statistical  analysis  of  tlie 
vertical  temperature  sounding  problem,  it  is  necessary  to  acquire  a valid 
"truth"  model  or  data  set  with  whicli  to  perform  comparison  studies.  The 
most  common  approach  is  to  use  radiosonde  dat.i,  wtiich  provides  relatively 
accurate  simultaneous  meastireraents  of  pressure,  liumidity  and  temperature 
and,  from  ground  tracking,  altitude.  Untortunately,  there  were  only  about 


50  radiosonde-DMSP  final  matches  using  the  3 hr.,  60  nm  criterion.  This 
number  is  much  too  small  for  performing  detailed  statistical  analysis. 

However,  these  data  can  be  used  to  infer  the  validity  of  the  TTPACK  data 
for  use  as  "truth". 

The  TTPACK  data  are  based  on  an  analysis  field  and  include 
numerical  interpolations  between  radiosondes,  including  the  latest  available 
measurements,  using  Cressman  weighting.  In  contrast,  the  DMSP  first 
guess  profiles  are  based  on  forecast  fields  and  do  not  include  the  latest 
available  data.  Another  advantage  of  using  TTPACK  data  is  that  it  uses  a 
model  which  partitions  the  ground  into  regions  (about  20  x 20  nm)  whereas 
the  radiosonde  data  are  point  measurements.  Since  the  DMSP  measurements 
are  over  regions  due  to  the  finite  resolving  power  of  the  sensors,  it  is 
more  appropriate  to  use  TTPACK  data,  if  it  assumed  that  the  TTPACK  errors 
relative  to  the  radiosonde  measurements  are  small. 

This  hypothesis  was  tested  using  the  data  supplied  by  AFCd.. 

The  results  are  summarized  in  Table  4.2.  Inspection  of  the  data  reveals 
that  the  TTPACK  profiles  match  the  radiosonde  profiles  more  accurately  than 
they  match  the  DMSP  final  profiles  and  are  more  accurate  than  the  match 
between  DMSP  final  and  radiosonde  profiles.  On  the  basis  of  statistics 
for  Feb.  14  to  Feb.  28,  latitude  45°N  -+  15°N,  TTPACK  is  an  unbiased  estimator 
at  most  levels  (perhaps  excluding  1000  mb,  30  mb  and  50  mb).  The  rms  per- 
formance is  'vl°K  at  levels  below  50  mb  witlt  the  exception  of  300  mb.  This 
is  consistent  with  the  fact  that  a RAOB  is  a point  measurement  while  TTPACK 
is  an  area  average.  In  addition,  RAOBs  are  generally  expected  to  have  a 

.5  to  1°K  rms  noise.  Exact  error  assignment  is  impossible  due  to  the  small 

s t 

number  of  samples.  This  should  be  contrasted  with  the  1--'  guess  - RAOB 
errors.  Large  mean  errors  exist  between  100  and  500  mb  and  1000  mb  (magnitude 
'\'1‘’,  often  2 or  .i^K).  This  trend  is  also  seen  in  the  DMSP-TTPACK 
errors.  The  sizes  of  these  samples  are  more  substantial.  These  means  should 
be  regarded  as  real.  The  same  general  trends  may  be  seen  in  the  data  for 
February  14  to  21,  latitude  62°N  to  90°N. 

A similar  trend  has  been  found  in  the  height  data  in  TTPACK-KaOB. 

The  mean  error  is  50  feet  below  150  mb.  Again  the  sample  is  small  (10). 

For  levels  below  250  mb  the  mean  error  is  less  than  10  feet.  Due  to  the 
fact  that  both  DMSP  and  TTPACK  represent  a spatial  average,  TTPACK  is  thus 
considered  to  be  a representative  ground  truth. 
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Data  matched  to 
within  4 hr,  90  rj- 


VALIDATION  OF  TTPACK  AS  "TRUTH"  DATA 


4.4  EXPERIMENTAL  RESULTS 

4.4.1  Statistical  Analysis  of  First  Guess  Temperatures 


and  Radiances 

The  rms  valuo.s  and  correlations  between  the  13  first  guess 
temperature  levels  and  the  8 IR  channels  are  given  in  Table  4.3  and 
the  associated  means  and  tlie  TTPACK  means  are  shown  in  Table  4.4.  Note 
that  channels  5 and  6 have  almost  identical  correlations,  implying 
that  they  are  almost  colinear  in  predictor  space.  This  suggests  that 
either  of  the  channels,  but  not  both,  should  be  used  in  the  retrievals. 

This  is  in  fact  presently  done  in  SOUNDER  by  using  only  channel  5. 

The  rms  first  guess  temperatures  vary  from  4.1  to  20.4°C  and 
tlie  rms  radiances  vary  from  3.1  to  24.9  ergs.  The  number  of  higli 
correlations  between  the  temperatures  and  between  the  radiances  is 
significant  and  suggests  that  there  is  a great  deal  of  redundancy  in 
each  type  of  predictor  data.  This  conclusion  is  verified  in  the  results 
presented  in  the  sequel.  Tlie  large  number  of  cross-correlations  between 
temperature  and  radiances  implies  that  there  is  a great  deal  of  redundancy 
between  the  two  types  of  predictor  data  as  well. 

As  discussed  earlier,  the  data  to  be  used  as  the  dependent 
or  "truth"  set  is  the  TTPACK  data.  Thus,  our  objective  is  to  reduce  the 
discrepancies  between  the  first  guess  and  TTPACK  profiles.  The  means  and 
covariances  for  tlie  first  guess  errors  (first  guess  - TTPACK)  are  given 
in  Tables  4.5  and  4.6  for  the  13  TTPACK  levels.  It  can  be  seen  that  the 
assumptions  used  in  the  minimum  information  solution  of  equal  diagonal 
elements  and  zero  off-diagonal  elements  is  subject  to  some  question.  Fleming 
and  Smith'  ‘ '’have  compared  the  performance  of  the  "full  statistical"  and 
the  minimum-information  methods  and  found  that  the  "full  statistical"  method 
was  superior.  They  attributed  its  superiority  to  the  fact  that  it  carries 
more  information  about  expected  statistical  variations  and  correlations 
in  the  predictor  data.  It  was  also  found  that  tlie  minimum-information 
accuracy  was  liighly  sensitive  to  the  accuracy  of  the  first  guess  profile. 

This  is  not  surprising  since  the  minimum-information  solution  will  generally 
be  (uniformly)  close  to  the  first  guess  profile. 


TABLE  4.4 


PREDICTOR  AND  TTPACK  MEANS-WEEK  #1 


Temperatures 

Radiances 

Level 

mb 

Mean 

First  Guess 

TTPACK 

Channel 

Mean 

1 

10 

-50.3 

o 

CO 

1 

1 

53.861 

3 

30 

-57.1 

-57.7 

2 

43.885 

4 

50 

-63,2 

-62 . 3 

3 

40.895 

6 

100 

-70.  3 

-71.9 

4 

51.601 

7 

150 

-63.4 

-62.5 

5 

68.514 

8 

200 

-55. 5 

-55,1 

6 

81.787  , 

9 

250 

-48.8 

-48.8 

7 

99.341 

10 

300 

-42.4 

-42.4 

8 

112.239 

11 

400 

-29.2 

-29.0 

12 

500 

-19.0 

-17.9 

13 

700 

-2.9 

-2.2 

14 

850 

4.3 

5.0 

1 

15 

1000 



10.9 

11.8 





1 

Temperatures  in 
Radiances  in  ergs. 


Mean  First  Guess  Errors 


MEANS  AND  COVARIANCES  OF  FIRST  GUESS  TEMPERATURE  ERRORS  - WEEK  #1 


First  Guess  Errors: 


MEANS  AND  COVARIANCES  OF  FIRST  GUESS  TEMPERATURE  ERRORS  - WEEK  #2 


r 


4.4.2  Discussion  of  Inversion  Methods 

A serie.s  of  experiment. s were  carried  out  to  evaluate  the  minimum- 
information  solution,  using  the  statistical  data  base  described  in  detail 
in  Section  A. 3.  The  statistical  analysis  consisted  of  the  following  three 
essential  steps: 

(i)  Develop  linear  regression  models  for  predicting  the 

temperature  profiles  from  the  first  guess  profiles  and  the 
measured  spectral  radiances.  This  provides  a performance 
baseline  for  linear  models.  This  includes  Important  trade- 
offs involving  the  number  and  type  of  predictor  (independent) 
variables  and  sensitivity  to  parameters,  such  as  latitude 
and  zenith  angle. 

(ii)  Develop  a hybrid  inversion  technique  in  which  the 

minimum-information  program  is  cascaded  with  a linear 
regression  model.  The  regression  model  uses  the  minimum- 
information  solution  as  part  of  the  predictor  set. 

(ili)  Analysis  and  comparison  of  techniques  to  determine  information 
handling  capabilities  of  tlie  minimum-information  method  and 
to  determine  possible  sources  of  modeling  error. 

The  basic  structure  and  Information  flow,  in  simplified  form, 
of  the  retrieval  methods  studied  is  given  in  Figures  4.2  and  4.3. 

Figure  4.2  sliows  the  information  flow  for  either  the  minimum-information 
or  regression  method.  The  basic  inputs  required  are  tlie  initial  forecast 
field  T^(p)  and  the  spectral  radiances  r(v).  The  output  T(p)  is  the 
estimated  temperature  profile.  The  performance  of  the  methods  is  evaluated 
by  a statistical  analysis  of  the  retrieval  error.  In  Figure  4.3  the  informa- 
tion flow  for  the  hybrid  method  employing  both  minimum  information  and 
regression  methods  is  shown.  The  inputs  to  the  regression  model  are  the 
spectral  radiances  r(v)  and  the  differences  between  the  initial  forecast 

field  T (p)  and  the  minimum-information  retrievals  T .(p).  As  with  the 
o mi 

other  methods,  the  performance  is  analyzed  by  statistical  analysis  of  the 
retrieval  errors. 

The  structure  and  information  flow  for  height  retrievals  is  similar 
to  that  for  temperature  retrievals. 
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Filjurc  4.1  Information  Flow  For  Minimum-Information  Tochniquc 


120 


The  hybrid  method  was  used  in  order  to  better  analyze  tlie  minimum- 
information  method.  If  tlie  minimum- inf ormation  method  was  a linear  processor, 
then  the  performance  of  the  hybrid  and  regression  methods  would  be  identical. 
Since  the  minimum  information  method  is,  in  fact,  a nonlinear  processor, 
performance  comparison  of  the  hybrid  and  regression  methods  allows  the 
efficiency  of  the  nonlinear  processing  to  be  measured,  at  least  in  a 
qualitative  way. 

The  coefficients  for  the  regression  and  hybrid  models  were  found 
using  the  stepwise  regression  method  outlined  in  Section  4.2.3,  using  the 
TTPACK  data  as  truth. 

If  the  set  of  predictor  (independent)  variables  is  denoted  by  the 
vector  X and  the  regression  coefficients  by  the  vector  a,  then  the  temperature 
estimate  T for  the  i—  level  is 
T 

T = a:x  + b. 
i 1 1 

with  b^  a constant.  Denoting  the  TTPACK  temperature  by  the 

regression  error  is 


= a^x  + b . 
1 1 


T 


Ti 


The  coefficients  a.,  b.  are  found  to  minimize  the  squared  errors  over  the 
11 

data  base  selected  for  study. 


The  temperature  levels  and  spectral  channels  used  in  the  DMSP  analy- 


sis are  given  in  Table  4.7.  Only  13  of  the  15  standard  pressure  levels  were 
used  in  the  analysis  since  TTPACK  for  20  mb  and  70  mb  were  not  available. 

The  window  channel  (835  cm  is  used  in  the  SOUNDER  program  to  provide  an 
estimate  of  surface  temperature  and  to  aid  in  the  cloud/no  cloud  decision 
process.  Tiie  ratio  of  spectral  radiance  in  channel  4 to  channel  7 is  used  in 
determining  whether  clouds  are  present  or  not.  Both  channel  7 and  channel  8 


radiances  were  used  in  the  regression  models. 


4.4.3  Statistical  Analysis  of  Inversion  Methods 

The  errors  in  the  first  guess  temperature  profiles  are  compared 


statistically  with  ttic  errors  in  tiie  minimum  information  retrievals  in  Table 
4.8.  The  mean  and  standard  deviations  of  errors  are  shown  for  both  week  #1 
and  week  If2  data.  I he  data  show  that  the  first  guess  is  biased  negatively 
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TABLE  4 . 7 


TEMPERATURE  LEVELS  AND  SPECTRAL 
CHANNELS  USED  IN  DMSP  ANALYSIS 


Temperature 

Level 

P (mb) 

CO, 

Channel 

u(cm  ^) 

Approximate  peak  of 

Transmission  Function  (mb) 

1 

10 

1 

668 

40 

2 

20  (1) 

2 

677 

60 

3 

30 

3 

695 

150 

4 

50 

(2) 

708'- 

350 

5 

70^') 

5 

725 

700 

6 

100 

6 

747 

1000 

7 

150 

2 

835^^^ 

8 

200 

8 

535(4) 

9 

250 

10 

300 

11 

400 

12 

500 

13 

700 

14 

850 

15 

1 

1 

1000 



Notes : 


(1)  not  used  in  analysis  (TTPACK  data  not  available). 

(2)  reference  frequency  used  in  minimum- inf ormation 
solution . 

(3)  window  channel 

(4)  water  vapor  channel 
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Level  I First  Guess  Final  ' First  Guess  Final 


COMPARISON  OF  DMSP  FIRST  GUESS  AND  FINAL  ERRORS 


at  nearly  all  levels  for  both  week  one  and  two.  The  DMfiP  final  errors  at 
the  same  levels  are  siqnif icantly  different:  the  bias  at  10  mb  stays  at 
about  the  same  maqnitude  but  reverses  sign,  while  the  biases  at  30,  50 
and  100  mb  are  all  negative  (between  -1.5  and  -4°C)  for  both  weeks. 

Since  the  minimum  information  iterations  terminate  on  the  basis  of 
radiance  matching  criteria,  this  implies  that  colder  than  actual  temper- 
atures are  required  to  match  the  measured  and  predicted  spectral  radiances. 
This  problem  could  bo  caused  by  errors  in  the  computed  transmittance 
functions  or  propagation  of  tropopause  errors.  Between  150  and  300  mb, 
the  performance  is  somewhat  different.  The  mean  first  guess  errors  are 
generally  less  than  one  degree;  however,  the  final  errors  are  between 
plus  2 and  4 degrees.  These  errors  could  be  due  to  compensation  of  the 
colder  than  actual  estimates  at  higher  altitudes  in  order  to  satisfy  the 
radiance  matching  threshold.  The  data  show  that  the  rms  error  in  the 
first  guess,  which  generally  varied  between  7 and  16°C,  was  significantly 
reduced  by  the  minimum  information  method  above  300  mb.  The  rms  DMSP 
final  errors  varied  bewteen  3 and  6°C.  Thus,  use  of  the  spectral  radiances 
significantly  reduced  the  variation  from  TTPACK  for  pressure  levels  less 
than  300  mb. 

For  pressure  levels  greater  than  300  mb,  the  situation  changes 
dramatically.  The  mean  first  guess  errors  are  relatively  small,  although 
they  appear  negatively  biased,  and  arc  generally  less  than  1°C.  The  mean 
DMSP  final  errors  are  essentially  the  same.  The  rms  errors  are  reduced 
generally  from  a range  of  3 to  8°C  to  a range  of  3 to  4°C. 
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In  order  to  test  the  statistical  significance  of  the  linear 
regression  models,  the  data  was  partitioned  into  two  sets.  Set  #1 
includes  all  soundings  from  February  lA  to  February  21,  1975.  Set  #2 
includes  all  soundinqs  from  February  22  to  February  28,  1975.  These 
two  data  sets  were  then  used  interchangeably  as  either  training  (dependent) 
or  test  (independent)  sets.  In  the  numerical  results  to  be  presented, 
the  notation  a ‘ b is  used  to  indicate  that  the  regression  model  is  derived 
from  set  //a,  wliile  the  regression  model  is  tested  on  set  #b.  All  four 
possible  combinations  (1*1,  1*2,  2*1,  2**  2)  are  utilized. 

The  RMS  regression  errors  obtained  by  using  all  13  first  guess 
levels  and  all  eight  spectral  radiances  are  shown  in  Figs.  4.4  and  4.7  for 
the  four  different  partitions  of  the  data  time  set  (1  *1,  1*2,  2 ■*  1 , 

2 ■>•2).  Shown  in  each  plot  is  the  effect  of  the  latitude  and  zenith  angle 
partitions.  The  "all  data"  curves  were  used  to  determine  a reasonable 
and  useful  bound  on  accuracy  of  temperature  retrievals.  That  is,  it  should 
be  possible  in  practice  to  realize  this  performance,  since  the  temperature 
retrievals  are  made  over  smoothed  data.  The  addition  of  structure  which  may 
be  added  via  partitioning  of  the  original  data  set  (into,  e.g.,  latitude  and 
zenith  angle  bands)  and  efficient  nonlinear  processing  would  be  expected 
to  yield  better  performance  since  the  variability  of  each  of  the  subsets 
is  less  than  the  variability  of  the  original  data  set.  It  is  interesting, 
as  will  be  subsequently  shown,  that  this  effect,  though  present,  is  not 
great  enough  to  invalidate  the  "all  data"  model  as  a possible  temperature 
inversion  method.  Several  conclusions  can  be  made  from  this  data. 

(1)  performance  is  almost  independent  of  which  training 

and  test  sets  were  used.  This  implies  that  both  weeks  have 
essentially  the  same  statistical  variability  and  that  regres- 
sion models  from  either  week  can  be  used  for  accurate  tempera- 
ture inversions  for  the  other  week. 

(2)  RMS  errors  for  all  cases  tended  to  be  higher  at  100  mb,  250  mb 
(near  tropopause)  and  at  1000  mb.  The  increase  at  100  mb  could 
be  explained  by  the  lack  of  a weighting  function  which  peaks 
nearby.  100  mb  is  approximately  midway  between  the  peaks 

at  60  mb  (for  the  695  cm  ^ cliannel)  and  150  mb  (for  the 
708  cm  ^ channel).  The  higher  errors  at  tropopause  are  not 


. >.  > 


Pressure  (mb) 


( 


unexpected.  Thie  relatively  high  errors  at  1000  mb 
are  probablv  related  to  the  errors  in  estimating  the 
surface  t imiperatiire. 

(3)  The  sensitivity  to  latitude  variations  is  much  greater 

than  the  sensitivity  to  zenith  angle  variations,  especi- 
ally at  high  altitudes.  Regression  errors  are  about  0.5°C 
below  50  mb  at  low  latitudes  and  rise  to  between  2.2°C  and 
3.2°C  at  the  higher  latitudes. 

The  regression  coefficients  for  all  week  //I  data  and  all  week 
#2  data  are  plotted  in  Figs.  4.8  and  4.9.  Continuous  plots  are  shown  rather 
than  just  the  points  in  order  to  better  display  the  trends.  The  coeffi- 
cients used  correspond  to  tlie  case  whore  one  first  guess  temperature,  at  the 
same  level,  and  all  8 spectral  radiances  were  used  in  the  independent  data 
set  for  estimating  each  Level. 

E.xamination  of  Figs.  4.8  and  4.9  reveals  several  things. 

Channels  2,3  and  7 contribute  most  to  tin.’  temperature  below  30  mb.  Note 
that  the  channel  1 coefficient  at  these  Levels  is  negative,  indicating  a 
compensation  for  the  skirts  of  the  transmission  functions  for  channels  2 
and  3.  The  first  guesses,  indicated  by  T^,  do  not  contribute  significantly 
to  these  estimates.  Channel  4 contributes  mostly  to  the  temperature 
estimates  between  200  and  300  mb.  Note  that  the  channel  4 coefficient 
between  50  and  150  mb  is  negative,  indicating  compensation  for  the  skirts 
of  the  tran.smission  functions  for  channel  3.  Note  that  the  initial  temperature 
is  significant  in  computing  the  100  mb  estimate.  This  may  be  obtained  by  the 
fact  that  none  of  the  transmission  functions  has  a peak  in  this  area,  as 
mentioned  previously.  Channel  5 contributes  mostly  to  the  region  between 
200  and  300  mb  and  also  compensates  for  the  channel  6 skirts.  Channel  6 
contributes  to  temperature  estimates  between  50  and  150  mb  and  compensates 
for  the  skirts  of  the  transmission  functions  for  channels  4,  5 and  7. 

Channel  7 contributes  mostly  in  the  region  from  200  to  400  mb  and  at  1000  mb, 
reflecting  its  usefulness  at  sensing  the  surface  temperature.  Channel  8, 
the  water  vapor  channel,  docs  not  strongly  contribute  to  any  level  but  has 
a weak  contribution  between  200  and  400  mb. 
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Figure  4.9  Regression  Coefficients; 
Week  2;  All  Data 


Ttio  K r.jtio.s  for  the  rocjro.ssioii  modei.s  cjivoii  in  4.8  and 

4.9  are  given  in  Tables  4.9  and  4.10.  Only  those  predictor  variables  which 
were  significant  at  the  5%  level  were  retained  in  the  regressions.  A zero 
value  indicates  tiiat  the  variable  was  not  significant  at  the  5%  level. 
Examination  of  the  data  reveals  several  things.  The  first  guess  tempera- 
tures are  relatively  insignificant  above  100  mb.  Between  100  mb  and  300  mb 
they  are  moderately  significant,  and  are  very  significant  between  400  mb 
and  850  mb.  They  are  seen  to  be  insignificant  at  1000  mb. 

Channel  1 is  significant  for  temperature  estimation  at  30  mb, 
which  is  close  to  the  weighting  function  peak.  Channel  2 is  not  consistent 
bu^  contributes  at  100,  150  and  200  mb  and  weakly  at  1000  mb.  Channel  3 is 
quite  significant  from  30  to  150  mb  and  surprisingly  significant  at  1000  mb. 
Channel  4 contributes  only  moderately  at  the  250  and  300  mb  levels.  Channel 
5 is  significant  at  150  mb,  mostly  due  to  its  role  of  compensating  for  the 
skirts  of  the  channel  3 transmittance  function.  Note  that  channel  5 also 
contributes  at  1000  mb.  Channel  6 contributes  consistently  at  50,  250,  300, 
700,  and  850  mb  levels  for  both  weeks.  Channel  7 contributes  significantly 
at  10  mb  but  only  marginally  at  1000  mb,  where  it  should  be  contributing  the 
most.  However,  Channel  7 contributes  very  significantly  in  the  interval 
between  250  and  400  mb.  Channel  8 (the  water  vapor  channel)  is  seen  to  be 
relatively  insignificant  for  use  in  temperature  estimation  at  any  level. 

A . A . 3 . ] i t ivi ^ to  Tra ipj-Pil  Test  Set  Data 

The  effect  of  using  different  training  and  test  sets  on  regression 
performance  was  evaluated  to  determine  the  statistical  significance  and 
predictive  power  of  the  regression  models.  The  results  are  summarized  in 
Eig.  4.10.  All  four  possible  combinations  of  training  and  test  sets  (1->^1, 

1^2,  2»1,  2 ^2)  wore  used.  When  the  training  and  tost  sets  wore  identica.l , 
the  rms  errors  were  generally  less  than  2°C.  However,  the  errors  at  100 
mb  and  250  mb  wore  slightly  larger,  about  2.5°C.  Tlie  largest  errors  occurred 
at  1000  mb,  and  wore  about  3.5°C.  Note  that  botti  curves  are  almost  identical, 
indicating  that  both  the  week  1 and  week  2 data  sets  have  the  same  statistical 
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TABLE  4.9 

F RATIOS  FOR  REGRESSION  MODEL  - V.TIEK 
• Use  only  one  first  guess  temperature  at  each  level. 
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TABLE  4.10 


F RATIOS  FOR  REGRESSION  MODEL  - WEEK  #2 
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• use  only  one  first  gue^s  temperature  at  each  level. 
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variability.  WlxtMi  diflorent  training  anti  test  sets  were  used,  the 
errors  were  slightly  degraded.  The  maximuin  degradation  in  rms  value 
is  0.b°C;  however  the  degradation  is  generally  less  than  0.4°C.  The 
shape  ot  these  curves  is  similar  to  that  ol  the  first  two.  The  conclu- 
sions to  be  drawn  from  these  results  are: 

(1)  Both  week  1 and  week  2 data  sots  are  large  enough  to 
be  statistically  significant.  That  is,  either  one  can 
form  tlie  basis  for  developing  a valid  regression  model. 

(2)  Either  of  the  regression  models  may  be  used  for  accurate 
temperature  retrievals  for  the  other  week. 

4.4. 3. 2 Comparison  of  Retrieval  Methods 

A scries  of  comprehensive  tests  were  made  on  the  minimum-informa- 
tion, regression  and  hybrid  methods  discussed  earlier  (cf.  Section  4.4.2). 
The  methods  tested  were: 

(a)  minimum-information. 

(b1  regression  method,  using  the  entire  first  guess  temperature 
profile  in  the  predictor  set. 

(c)  regression  method,  using  only  the  first  guess  temperature 
at  the  level  of  interest  in  the  predictor  set. 

(d)  hybrid  method,  using  the  entire  first  guess  temperature 
profile  in  the  predictor  set. 

The  results  of  these  tests  are  summarized  in  Figs.  A. 11-4. 32,  in 
which  the  retrieval  errors  in  “C  are  plotted  for  all  data  and  the  two 
latitude  and  zenith  angle  partitions.  The  most  significant  results  may  be 
summarized  as  follows: 

(1)  The  minimum-information  retrievals  are  uniformly  worse  than 
the  other  three  methods  when  the  regression  and  tests  are  Identical,  with 
the  exception  of  the  2 -* 2 case  for  y > 45“  at  250  mb. 

(2)  More  s i gn  i f leant  1 v,  tiie  minimum-information  retrievals  are 
uniformly  worse  than  the  other  methods  when  different  training  and  test 
sets  are  used,  with  the  exception  of  the  case  y >45°  near  300  mb.  The 
performance  degradation  is  quite  severe  in  many  cases,  especially  above 
100  mb,  where  it  ranges  generally  from  2'^C.  to  6°0.  The  degradation  for 
pressures  below  100  mb  are  generally  about  l°l'. 
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(3)  Tlic  fall  regression  raetiiod,  using  tlie  entire  first 

* 

guess  profile  (all  3'^ ) . is  uniformly  better  tlian  the  other  methods. 

This  result  is  not  unexpected  when  tlie  training  and  test  sets  are  identical. 
However,  tliis  result  holds  also  for  tiie  casi’S  when  the  training  and 
test  sets  are  different.  The  rms  errors  were  generally  about  1“C  above 
50  mb,  between  2 aitd  3°r  at  100  mb  and  250  mb  and  between  1 and  2°t:  below 
300  mb.  Errors  were,  however,  somewhat  larger  at  1000  mb,  as  higli  as  5°i.'. 

This  is  annarently  due  to  tlie  wav  in  which  t'iie  dear-column  measurements 
are  nicked.  .Measurements  may  be  made  even  when  cloud  cover  is  present. 

(4)  Tlie  effect  of  eliminating  all  but  one  first  guess  tempera- 
ture level  has  little  ef feet,  generally  0.2°r.  or  less.  There  were  several 
cases  for  which  accuraev  was  significantly  degraded  - up  to  1.5“(’  at  100  mb 
(2  ->•  2,  all  data)  . 

(5)  Perhaps  the  most  surprising  result  is  that  the  performance 

of  the  hybrid  method,  denoted  "minimum-information  residuals"  was  uniformly 
wo_^se  than  the  regression  method.  The  implication  of  this  is  that  the 
minimum-information  method  is  nonlinearly  processing  the  data  in  a non- 
optimum  manner.  It  appears  that  information  is  somehow  lost  in  the  minimum 
information  solution.  As  discussed  previously,  if  the  nonlinear  processing 
was  efficient,  the  hybrid  method  would  be  expected  to  achieve  the  best 
performance.  The  degradation  of  the  hybrid  method  is  quite  striking  in  some 
cases  - often  as  high  as  1 '’C  and  even  reaching  2°P. 

It  should  be  noted  that  the  rms  errors  for  the  hybrid  solution 
are  quite  sensitive  to  zenith  angle  above  200  mb,  with  significant 
degradation  from  the  optimal  (regression)  results  appearing  at  large  zenith 
angles,  but  little  degradation  at  low  zenith  angles.  Since  the  zenith 
angle  computations  arc  significant  only  for  high  values,  this  implies  that 
the  zenith  angle  processing  in  the  minimum-information  method  is  degrading 
the  content  of  the  data.  Further,  since  the  hybrid  and  minimum-information 
performance  is  quite  similar  in  this  region  at  G>0.7,  minimum- information 
appears  to  be  doing  a good  job  of  linear  data  processing.  Thus,  the 
implication  is  that  the  information  is  being  lost  in  the  nonlinear  processing. 
It  is  possible  that  program  modifications  to  eliminate  the  zenith  angle  pro- 
cessing would  enhance  the  minimum-information  solution  performance. 

*Slight  exceptions  such  as  for  the  case  y >45°,  1 >■  1 at  150  mb  can  be 
attributed  to  the  stepwise  regression  method  and  the  way  variables  are 
selected . 
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(f))  Thf  minimum-information  rms  orrors  are  significantly  larger, 
iij)  to  f)°K,  ubovi’  100  mh  at  high  Jatitucies  than  at  low  latitudes.  This  may 
be  .itfribiited  to  the  larger  variation  in  temperature  profiles  at  high  lati- 
tudes, the  increased  v.iriability  in  the  meastirements  dvie  to  cloud  presence 
and  the  method  in  wliich  the  first  guess  profile  is  computed.  i'i)r  pressure  level 
alx've  10  mb  ami  latitudes  ^ 60°N  only  two  guess  profiles  are  used,  one 
for  summer  and  one  for  winter. 

(7)  The  minimum-information  errors  are  relatively  independent 
of  zenitii  angle,  althougli  they  are  somewhat  larger  between  50  mb  and  200  mb 
for  0 < 0. 7 rad  ian . 

(0)  The  reqression  rms  errors  vary  somewhat  with  latitude. 

Below  45°N,  the  errors  atwvo  50  mb  are  small  (about  0.5  K)  but  increase 
to  Ix'tween  2°  ar.d  3“k  at  100  mb.  Al>>ve  45''N,  the  errors  are  generally 
hetwi-en  2°  and  at  all  levels.  The  increased  error  ar  100  mb  at  Low 

latitudes  is  probably  due  to  the  cjuiescent  nature  of  the  temperature 
profiles,  which  makes  it  difficult  to  obtain  estimates  in  the  middle  of 
transmittance  function  peaks. 

(9)  Tlie  regression  errors  are  insensitive  to  changes  in 
zenith  angle.  The  results  obtained  using:  (1)  all  data,  (2)  0 < 0.7 
and  (3)  G > 0.7  are  almost  identical  in  all  four  cases.  Since  the  zenith 
angle  was  not  included  in  the  regression  model,  this  implies  an  essential 
Independence  of  the  information  with  respect  to  zenith  angle.  Thus  tlie 
zenith  angle  does  not  appear  to  carry  information  relevant  to  the  problem 
of  temperature  sounding.  This  also  implies  that  the  zenith  angle  computations 
in  SOUNDER  may  not  be  required,  althougli  a more  careful  study  of  the  inter- 
action between  0 and  other  internal  program  variables  would  have  to  be  made 
to  corroborate  this  conclusion. 


Pressure  (mb) 


Temp)erature  Inversion  Error  Comparisons 


01  23456  78  9 10 

rms  Error  (°K) 


Figure  4.11  Temperature  Inversion  Error 
Comparisons;  All  data(l*-l) 
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Figure  4.16  Temperature  Inversion  Error 
Comparisons;  All  Data  (l-*-2) 
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Figure  4.17  Temperature  Inversion  Error 
Comparisons;  X < 45°  (1>2) 
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Figure  4.18  Temperature  Inversion  Error 
Comparisons;  A > 45°  (l->2) 
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Figure  4.20  Temperature  Inversion  Error 
Comparisons;  0 > 0.7  (l->-2) 
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Figure  4.21  Temperature  Inversion  Error 
Comparisons;  All  Data  (2'-l) 
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Figure  4.22  Temperature  Inversion  Error 
Comparisons;  X < 45*  (2-vl) 
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Figure  4.30  Temperature  Inversion  Error 
Comparisons;  6 > 0.7  (2->2) 
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Figure  4.J2  Comparison  of  Mean  Tempera- 
ture Inversion  Errors;  All  Data;  Week  2 
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A . A . 3 . 3 Pressure- Inver s ion 

rnvcrsiens  wore  porformod  using  prossiiri'-he  ight  data  for 
week  i/ 1 . The  data  was  idontic'al  to  tiial  used  for  t onipera  t ure  inversions. 

The  results  of  perfoiming  stepwise  regressions  on  the  data,  using  all  first 
guess  temperatures  and  speetral  radiances  in  t lie  predictor  set,  are  shown  in 
Table  4.11. 

A . A . 3 . A Comparison  of  Temperature . Brightness  Temperature 
and  Height  Inversions 

Linear  regression  models  were  also  generated  for  brightness 
temperatures  and  for  pressure  heiglit  data.  The  brightne.ss  temperature 
inversions  were  obtained  by  first  converting  the  first  guess  and  TTPACK 
profiles  to  brigiitness  temperature  via  Planck's  equation  and  then  developing 
models  to  minimize  tlie  mean  squared  errors  between  them,  using  the  brightness 
temperatures  and  spectral  radiances  as  predictors.  The  retrieved  brightness 
temperatures  were  then  .inalyzed.  The  pressure-lie  ight  data  used  corresponded 
exactly  to  the  cases  used  tor  temperature  inversions.  An  intercomparison 
of  the  three  types  of  retrieval  errors  was  made  using  the  percentage  of 
total  variation  explained,  adjusted  for  the  number  of  degrees  of  freedom. 

The  results  are  shown  in  Table  4.12.  The  three  types  of  predictands  shown 
are 

(1)  temperature  (“O 

(2)  brightness  temperature  (ergs),  obtained  by  converting 

temperatures  to  equivalent  radiances  via 
Planck's  equation 

(3)  heights  (ft) 

Ttu'  results  tor  the  temperatures  and  briqhtness  temperatures  arc  quite 

similar  and  show  a marked  deqiadation  around  tropopause  (200  mb).  The 

tieinlit  predictands  show-  no  deqradation  at  this  level  since  there  is  no 

sharp  discontinuity  in  sbipt'  as  is  the  ca.se  for  temperature  at  tropopause. 

The  heiqht  predictands  explain  alout  dOf  of  the  variance  Ijotween  100  mb 

and  500  mb,  while  the  temperature  predictand.s  qenerally  explain  between 

HU.  and  dPA  in  the  same  reqion,  except  lU'ar  troi'opauso  where  the  value 

2 

cirop.s  to  alx-iiit  70".  Note,  however,  that  R for  heiqhts  drop.s  liolow  the 
temperature  values  below  700  mb.  This  is  p'rolxibly  due  to  the  decreased 
uncertainty  in  heiqht,  or  thickness,  first  quess  profiles  near  the  earth. 
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TABLE  4.11 


REGRESSION  ANALYSIS  FOR  HEIGHTS  - WEEK  #1 


TABLE  4.12 


2 

Percentage  of  Variation  Explained  (R  ) for  Temperature, 
Brightness  Temperature  and  Height  Inversions 
All  Data  - Week  1 
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4 . 4 . A Regression  Coe f f icients 

The  regression  coefficients  for  both  week  1 and  week  2 models  are 
given  here.  Tlie  models  selected  for  display  utilize  only  the  first  guess 
temperature  at  the  level  being  retrieved  (the  1 case),  since  these  gave 

essentially  the  same  performance  as  using  the  full  profile.  Thus  the  predictor 
set,  before  stepwise  regression  contained  a single  first  guess  temperature  (in 
degrees  C)  and  eight  spectral  radiances  (in  ergs).  The  regression  equation  for 
the  1—  level  is 


T,  = a. .T  . 
1 11  oi 


8 

+ I 

k=l 


ik  k 1 


i = 1.2, 13 

2-1 

with  T.,  T Y.  measured  in  °C  and  r,  measured  in  ergs/ (cm  ‘sec*  ster 'em  ). 

1 oi  1 k 

The  coefficients  obtained  using  all  data  are  displayed  in  Tables  4.13  and  4.14 
for  all  13  temperature  levels.  It  is  interesting  that  the  regression  models 
are  somewhat  different  tor  the  two  weeks  although  the  same  general  trends  are 
present,  as  would  be  expected.  The  fact  that  both  models  yield  quite  accurate 
retrievals  for  either  week  indicates  that  further  reduction  in  the  number  of 
predictors  is  probably  possible  with  very  little  loss  of  retrieval  accuracy. 
Further  work  needs  to  be  done  in  the  area  of  model  reduction.  Possible 
methods  which  should  be  explored  Include  the  statistical  eigenvector  method  and 
the  use  of  coefficient  shrinking  techniques  such  as  ridge  regression  and  Stein's 
method  (see,  e.g.  Dempster,  et  al. 
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1000  i ! .368  I 0 .668  -1.413  | 0 .533  0 .250  0 -25. 


REGRESSION  COEFFICIENTS  FOR  WEEK  ?/ 


.867  -1.423  ' 0 1.373  -.857  ' .689  0 -60.6 


5 CONCLUSIONS 


■ 

1.  On  the  average,  radiances  measured  at  all  latitudes  and  longitudes 
for  Bird  8531  correspond  closely  with  those  from  Bird  9532. 

2.  Radiances  measured  at  all  Bird  8531  channels,  except  4,  are 

greater  (warmer)  than  those  from  Bird  9532.  * 

3.  Criteria  used  to  insure  clear  column  viewing  by  discriminating 

against  low  clouds  and  hot  surfaces  (CLDDIF  and  HOTSF)  are  far  ^ 

more  restrictive  in  today's  operations  at  AFGWC.  Results  show 

that  the  present  SOUNDER  accepts  only  1/4  as  many  cases  for 

completing  an  inferred  profile  as  it  did  in  the  197  5 research  data 

set.  In  fact,  out  of  10  cases  with  nearly  coincident  satellite  and 

RAOB  soundings  only  one  case  could  be  retrieved  using  present- 

day  criteria.  Many  of  the  most  significant  meteorological 

tem.perature  profiles  are  being  eliminated. 

4.  Signal  to  noise  ratio  and  the  epsilon  value  (SNOIS  and  EPS)  are 
used  to  control  convergence  for  a solution  and  to  mionitor  the 
maximum  difference  allowable  between  calculated  and  measured 
radiance  from  any  channel.  Wide  variations  in  SNOIS  and  EPS 
produced  insignificant  differences  in  the  retrieved  temperature 
profile.  Tne  additional  computei  time  required  to  obtain  convergence 
of  a solution  was  not  warranted  by  the  use  of  the  more  restrictive 
values  of  SNOIS  and  EPS. 
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5.  Variatious  in  the  water  vapor  profile  produce  insignificant  effects 
on  retrieved  temperatures  in  high  latitudes  but  become  dominant 
for  mid-latitude  and  tropical  regions.  Errors  due  to  water  vapor 
effects  arc  most  pronounced  on  850  and  700  mb  retrieved  tempera- 
tures . 

6.  Underlying  surface  effects  influence  the  measured  radiances  at  all 
channels,  even  tnosc  whose  weighting  functions  are  near  zero  at 
the  surface.  Accounting  for  the  underlying  surface  effects  is  more 
important  than  water  vapor  effects  in  the  mid-and  high-latitudes . 

The  type  and  intensity  of  the  surface  effects  are  highly  dependent 
upon  the  low  level  air  temperature  profile.  In  fact,  a good  correla- 
tion (0.9  12)  exists  between  the  low  level  air  temperature  gradient 
derived  from  the  forecast  data  and  the  retrieved  temperature  errors 
at  1000  mb  level.  It  appears  that  significant  improvements  in 
retrieval  of  1000  mb  air  temperatures  are  attainable. 

7.  Detailed  statistical  analysis  of  DMSP  minimum-information  re- 
trievals has  been  made  using  AFGL  data  irom  February  14-28,  197  5, 
with  TTFACK  data  used  as  truth. 

8.  A series  of  linear  regression  models  were  developed  using  forward 
stepwise  regression  on  TTPack  data  as  the  truth  model.  First  guess 
temperatures  and  spectral  radiances  from  eight  DMSP  SSE  channels 
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were  used  in  the  predictor  set. 


Regression  metliod  yields  significantly  lower  mean  and  rnis  errors 
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9 . 

tnan  the  minimum  information  technique  over  the  entire  Northern 
Hemisphere  and  at  all  zenith  angles. 

10.  Minimum-information  method  appears  to  degrade  information  during 
nonlinear  processing  of  input  data,  (e.g.,  temperature  dependence 
of  water  vapor  transmittance,  CO;  transmittance  function  aependence 
on  zenith  angle). 

11.  Zenitii  angle  does  not  afipeai  to  he  a significant  predictor. 

12.  It  appears  that  simple  rcgres.aon  models  may  be  fitted  on  weekly 
data  and  used  for  accurate  temperature  retrievals  for  the  following 
week.  Tins  implies  that  the  data  base  contained  enough  statistical 
variation  to  yield  statistically  significant  results.  More  importantly, 
the  results  indicate  that  the  regression  models  might  be  used 
operationally  for  temperature  retrievals. 
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6 RECOMMENDATIONS 


1.  The  acceptance  criteria  to  distinquish  between  cloudy  or  hot  surfaces, 
CLDDIF  and  HOTSEC,  need  to  be  reviewed  because  underlying  surface 
effects  appear  to  bo  causing  vertical  soundings  to  be  discarded  for 
low  nadir  angles  and  for  the  most  significant  types  of  meteorological 
profiles . 

1.  Underlying  surface  effects  must  be  defined  and  incorporated  before 
significant  analytical  or  statistical  operational  improvements  can  be 
expected  in  any  satellite  temperature  retrieval  near  the  earth's  surface. 

3.  A study  should  be  performed  using  the  ALLRADIANCE  data  to  investi- 
gate cloud  effects  and  reduce  their  effect  on  temperature  inversion 
errors . 

4.  The  effects  of  sensor  center  frequency  shifts,  water  vapor  uncertain- 
ties, and  transmission  function  variations  needs  to  be  thoroughly 
studied . 

5.  GAC's  Research  version  of  SOUNDER  is  an  effective  tool  that  should 
be  used  to  define  more  precisely  the  effects  on  and  improvements  to 
the  AFGWe  operational  version  for  retrieving  temperature  profiles  from 
DMSP. 

6.  Use  of  statistically-based  models  should  be  investigated  further  for 
other  seasons  and  with  a larger  data  base. 
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7.  'Hie  otfcct  of  zenith  angle  on  temperature  inversion  should  be 
thorougiily  investigated . 

8.  Based  on  this  study,  statistical  models  slrould  be  strongly  con- 
sidered for  implementation  at  AFGVVC. 

9.  The  use  of  additional  data,  such  as  microwave,  and  sea  surface 
temperatures  should  be  considered  for  the  temperature  inversion 
problem . 
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